Asymptotic Morse theory for the equation Av = 2v x A v y 
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abstract. - Given a smooth bounded domain f2 C R 2 , we consider the equation Av = 2v x A v y in il, 
where v : fl — > R 3 . We prescribe Dirichlet boundary datum, and consider the case in which this datum 
converges to zero. An asymptotic study of the corresponding Euler functional is performed, analyzing 
multiple-bubbling phenomena. This allows us to settle a particular case of a question raised by H. Brezis 
and J.M. Coron in 
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1 Introduction 

Let fl C R 2 be a smooth bounded domain. We shall denote by v,g two maps such that v : — > ]R 3 and 
g : dfl — > M 3 , with g smooth. Consider the problem 



(1) 



Av — H(£, v, Vu)w x A v y in f2, 
v = g on dft, 



where H is a smooth scalar function, v x , v y are the x and y-derivatives of v, £ = (x,y) and A denotes 
the cross-product in R 3 . 

Equation ([!]) has been the subject of several works, see for example the survey paper by K. Steffen 
and the recent paper Jl0| ]. Existence of solutions of (Q) when g = strongly depends on the topology of 



the domain. In fact we show using a Pohozahev-type identity, see Proposition 3.1, that equation (|lj) has 
no solution in any simply connected domain when g = 0. When i?(£,t>, Vt>) = H, a non-zero constant, 
such a result was proved by H. Wente, llfj], using reflection techniques and the Kelvin transformation. 
In the same paper, Wente also showed that if f2 is an annulus then the study of (0) can be reduced to an 
ordinary differential equation and (|l|) does have a non-trivial solution when v = on 90. Thus equation 
(Q) presents features similar to the Yamabe equation on domains with Dirichlet boundary condition, 
studied in particular by A. Bahri and J.M. Coron, 0]. In fact, part of the difficulty in studying (|l|) is 
that it is invariant under conformal transformations. This invariance forces the associated variational 
problem to exhibit non-compactness phenomena, like in the Yamabe problem on domains. We point out 
that in our case, contrary to the Yamabe problem, simply connected domains always admit only trivial 
solutions. For the Yamabe problem in dimension greater or equal than three, there are indeed examples 
of contractible domains which admit non-trivial positive solutions, see psj . 

From now on we consider the case of constant H, precisely v, V«) = 2. So problem ([l]) reduces 

to 



(2) 



Av — 2v x /\v y , in fi, 
v = g on 9f2. 



1 E-mail addresses: chanillo@math.rutgers.edu (S. Chanillo), malchiod@ias.edu (A. Malchiodi) 



Under the assumption ||g||oo < 1, S. Hildebrandt, |jl8| , constructed a solution of (Q) with minimal energy 
called the small solution, while Brezis and Coron, K.Steffen, and M. Struwe, [|9|, constructed a 
second solution, referred to as the large solution. We remark that the assumption \\g\\oo < 1 is sharp, see 

0- 

Results similar to those regarding the Dirichlet problem hold for the Plateau problem, in which one 
looks for solutions of Au — Hu x A u y which are conformal and which map the boundary to a given curve 
(with free parametrization). As a result one obtains surfaces with constant mean curvature. 

We mainly focus on the following problem 



(3) 



Av — 2v x /\v y , in f2, 
u = eg on dil. 



We will study (^) turning it into a variational problem. In view of the non-existence result in |L7|| , it 
is natural to assume that the boundary datum is small. T. Isobe in particular, p(|-p2|, analyzed the 
behavior of the large solutions of Brezis and Coron in the limit e — > (the small solutions converges to 
the trivial one v = as e — ► 0). 

Let g denote the harmonic extension of g in Q, i.e. 



Ag = in Q; 
g = g on dCl. 



(4) 

If v is a solution of (^) and if we set v = u + e g, the function u solves 

\Au = Av = 2{u x +e g x ) f\{u y +e g y ) in fi; 
1 u = on dfl. 

Problem (^|) admits the Euler functional I e : -ffg(f2;R 3 ) — > M, which has the following expression 
(5) I e {u) = - / |Vu| 2 + \ { u ■ (u x Au y ) + e [ u- (u x A g y + g x A u y ) + 2e 2 [ u ■ (g x A g y ). 



J') 

The aim of this paper is to develop a Morse theory for the functional I s when e is small. In order to do 
this we take advantage of the perturbative approach in [jlj . We first recall from |^] that the fundamental 
solution {bubble) of the equation 

(6) Au = 2u x /\u y , inM 2 

is the stereographic projection ir : M 2 — > S 2 C R 3 

2x 2y x 2 + y 2 - 1 

1 + x 2 + y 2 ' 1 + x 2 + y 2 ' 1 + x 2 + y 



Our analysis will use translations, dilations and rotations of the function in (Q) and we set 

(8) RS a ,x(x, y) = Ro tt{\{x - Oi, y - a 2 )), 

for R £ 50(3), a = (01,02) £ M 2 and A > 0. The functions R8 a ^\ are mountain-pass critical points of 
the functional 

(9) I( u ) = l [ \S7u\ 2 + l [ u-(u x Au v ), u£V, 

1 JR 2 3 J R 2 

where T> denotes the functional space 

V=[u£Ll r (M 2 ;« 3 ) ■ \\"\\h= i IV.f - / j" 1 , ■ x 
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The space T> coincides with H 1 (S 2 ;M 3 ) after inverse stereographic projection. We point out that the 
functionals I and I £ are well defined and smooth on T> and Hq (n, R 3 ) respectively, see Section ^. It turns 
out that the manifold constituted by the <5's is non-degenerate for the functional / (modulo constants), 
as proved in pl[ Lemma 5.5, using an isoperimetric inequality. Proving the non-degeneracy condition is 
equivalent to classify the solutions of 

(10) Aw = 2 (w x A 5 y + S x A w v ) in R 2 , 

which is the linearization of (||) around 5, and to show that the only solutions are the tangent vectors to 
Z at 8, see equation (|l6|). We remark that equation (^) admits solutions of the form ir(z k ) (in complex 
notation) for any integer fc, see We will refer to these solutions as higher degree bubbles. For this 
reason we give in the Appendix an alternative proof of the non-degeneracy, which we believe could adapt 
naturally to the higher-degree case. 

To analyze the problem in f2, it is convenient to consider the functions PS = 5 — ip, where (p is defined 
in (pl|). PS is the element of i?Q(f2;R 3 ) closest to S in the Dirichlct norm. We may write 

k 

(11) 7i = ^Pi? l( 5 p! , Al +w, 

where Ri £ 0(3), A^ > 0, pi <E R 2 for all i, and w is orthogonal to the manifold J3 i=1 PRiS Pu Xi- Once 
we have the non-degeneracy property for J, then it is standard to prove that for suitable values of a and 
A also the manifold of projected bubbles is non-degenerate for I £ , and the same holds true for a finite 
sum of bubbles. This property allows us to solve the equation I' E (u) = in w (see Proposition |4.3| ), and 
thus our problem is reduced to a finite-dimensional one which involves an auxiliary functional I e (z) (see 
Section ^ depending only on {pi}i,{h}i and {Ri}i. Substituting (11) into I e and letting e — ■> 0, we 
expand I £ (z) for large values of Xi (roughly of order e _1 ). 

The large solution of Brezis and Coron corresponds to a one bubble solution when e — > 0, and has 
been studied in detail by T. Isobe, [E0)-||22]. However, from Theorem 0.3 in M it is clear that a more 
complicated configuration may occur. Thus to manufacture this type of solutions we are naturally led 
to a variational analysis of the functional (p) for multiple bubbles. We point out that from the work of 
Brezis- Coron the bubbles will not necessarily be all of degree 1. However the variational analysis is more 
difficult if we allow bubbles of arbitrary degree, and we will return to this point in a subsequent article. 

To state our results we need some notation. Given (a,£) = ((ffli, 0,2)1 i x i v)) G ^ x J7, let h\,h,2 '■ 
51 x — > R be the solutions of the problems 



(12) 



A^i(a,£)=0 inO, J A e /i 2 (a,£) = in n, 

>( a >6 = |=^ ondtt; \h 2 (a,C) = fff^l on d£l, 



see Remark 5^ (6). If G(a, £) denotes the Green's function of Q,, normalized so that G(a, £) ~ — log |a — £| 



for a ~ £, and if H(a,£) denotes the regular part of G (G(a,^) = — log |a — — H(a,t;)), then we have 
(13) /ii(a,£) = — ; /i 2 a,£) = — 5 . 



Let also 

(14) H(a) 



dhi dh 2 
dx dy 



It has been proved in [M-|22 (see also [^6|) that the function H plays a crucial role in studying the 
location of blowing-up solutions of (||), when the boundary datum converges to zero. In fact, the function 
H appears in t he e xpansion of I e (u), when u is of the form ( |ll| ) with fc = 1, as a self-interaction term, 
see Proposition |5.l| . The expansion for fc = 1 is essentially performed in the works of Isobe, but we derive 
it in a framework which is convenient to treat the case of fc > 1, see Section 0. We have the following 
result, regarding the function H. 
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Theorem 1.1 (a) For a simply connected domain Q there holds 

H(a) = 2e 2H{a a \ a e fi, 

where H(a,£) is the regular part of the Green's function o/O. 

(b) For a multiply connected domain Q there not exist in general a function F such that H(a) — 
F(H(a,a)). In particular, for some annulus of the form {p^ 1 < \z\ < p} , p > 1, the critical points of 
H(a,a) and of H(a) do not coincide. 



Theorem LI is proved in Section ||. The function H{a, a) is called the Robin function of the domain fl, 



see |6[. In dimension 2 it also appears in extremal problems related to the Moser-Trudinger inequality, 
where the critical points of H(a, a) are shown to be related to the conformal incenter of fl. Isobe showed 



that H > on any domain, see Remark 5.3, but did not analyze it further. Since H is defined by means 



of second derivatives of H, we need to use a global argumen t (th e Riemann mapping theorem) to compare 



the two functions H and H. The new feature of Theorem 1.1 is that the Robin function plays a role in 
concentration phenomena only for the case of simply connected domains. 

The regular part of the Green's function plays an important role in many problems with critical 
exponent in dimension larger than two, see ^j, ||, 0, [j^, p^| , [p5[ . The difference here is that the 
regular part does not appear directly in the expansion, we find the above function H instead, and we 
recover the regular part from the Riemann mapping Theorem. The Robin function is also related to the 
notion of conformal incenter, sec fllT 



The expansion of I £ (u) for multiple bubbles is performed in Section ^, see Proposition 7.4. It turns 



out that when Xi ~ e for all i, the mutual interaction among the bubbles is of the same order as the 
interaction with the boundary (through both the geometry of fl and the datum g). We observe that the 
interaction among the bubbles depend on their mutual orientation. 

There is a by-product of the expansion in Propositio n \l.4\ . It allows us to settle a particular case of 
a question raised by Brezis and Coron, see Section |^. In || the authors consider a sequence of solutions 
u„ of (|l]) and a sequence g n of boundary data which converge to zero in H^{dfl) n L°°(dfl). Under 
these conditions they prove that the sequence u n splits into a finite number of bubbles, and their image 
converge to a finite and connected union of spheres of radius 1. They ask whether every configuration of 
spheres can be obtained as a limit of solutions u n for a suitable sequence of boundary data g n . We have 
an affirmative answer if all the spheres pass through the origin. 

Theorem 1.2 Let D denote the unit disk in R 2 ; and let A = {Si U • • ■ U Sk} be any configuration of unit 
spheres, each passing through the origin o/R 3 . Then there exist a sequence g n : S 1 — > R 3 and a sequence 
of functions u n solving 

{Av = 2v x A v y in D, 
v — g n on dD, 

such that the image of the function u n converge to A in the Hausdorjf sense. 



The functions u n in Theorem 1.2 are constructed studying the interactions of the bubbles (of degree 1) 
with the boundary datum and among themselves. Choosing boundary data with an appropriate strong 
concentration at k points on dD, we show that the self interaction among the bubbles becomes negligible. 
Hence we can find solutions u n which are highly concentrated at k points close to the boundary of D and 
with prescribed orientations in R 3 . We remark that the order of concentration, roughly the parameter A 
in (||), turns out to be the same for all the bubbles. 

The case of spheres not passing through the origin is not treated here. We believe that it could 
be possible to achieve such configurations by considering bubbles with higher degree. In fact, in the 
recent paper by A. Bahri and S. Chanillo, ||, the authors showed that when considering changing-sign 
solutions of the Yamabe problem, the bubbles can exhibit different orders of concentration. If there is 
an analogy between higher-degree bubbles and changing-sign solutions of the Yamabe equation, then 
one could obtain bubbles with higher and higher concentration and with image not passing through the 
origin. This will be the object of a future work. 
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2 Notation and preliminary facts 

In this Section we introduce some notation and preliminary facts in order to tackle problem 
In the following D will denote the unit disk in R 2 

D = {£=(x,y)eR 2 : x 2 + y 2 < 1} . 

Let / : V — > R be defined by From |T^] the last term in I(u) is well defined on V, together with its 
Frechet derivatives. This makes I a smooth functional on T>. The same argument provides regularity of 
the functional I e on i?g(0;R 3 ). 

Using a finite-dimensional reduction, we are going to treat the functional I £ as a perturbation of /. 
In order to do this, it is essential to consider the critical points of the functional /, namely the solutions 
of 

(15) Au = 2u x Au y inR 2 , u G V. 

The stereographic projection ([?]) is indeed a solution of (|l5|), which we call fundamental solution or bubble. 
By invariance its translations, dilations and rotations are also solutions of ([IB]). We set 

(16) ~Z = {R5 a ,x{-) = RSx(--a) : A > 0, a € R 2 , R e SO(3)} . 

We remark that, since SO (3) is a three-dimensional manifold, Z is a six-dimensional manifold. 
We list now some useful expressions. Note that the function S (A (£ — a)) has the explicit form 



*(A(£-o)) = 



2A(f-a) A 2 |<e — a| 2 - 1 



,l + A 2 |e-a| 2 ' A 2 |e-a| 2 + l 
from which, if a is bounded away from 9f2, one can deduce 

(17) ^(A(^-a))~0|^,l-^| ? ^)+O(A- 3 ), for A large. 

Writing for brevity 5 instead of <5 a ,A> we compute some derivatives of 5. We emphasize that throughout 
the paper, unless explicitly stated, the point a will always be bounded away from dfl, namely we will 
assume dist(a, dCl) > to for some fixed To > 0. We have 

l + X 2 (y 2 -x 2 ) _ 4\ 3 xy _ 4\ 2 x 



'{l + X 2 (x 2 +y 2 )) 2, v XJZ (1 + A 2 (a; 2 + y 2 )) 2 ' v XJ (1 + A 2 (a; 2 + y 2 )) 2 ' 
(18) 

4A 3 xy 1 + AV-y 2 ) 4A 2 y 

[ y,1 ~ (l + \ 2 (x 2 +y 2 )) 2 ' [y,2 ~ (l + A 2 (x 2 +y 2 )) 2 ' 1 a j3 ~ (1 + A 2 (a; 2 + y 2 )) 2 ' 

From the last two formulas we deduce 

, . r . 8X 3 x , 8X 3 y 

( 19 ) fe A ^ = - (1 + ^ + ^ i ^ A ^ )2 = - (l + A 2 (x 2 +y 2 ))3 ; 



5 



The functions R8 at \\n do not belong to iJQ(f2;R 3 ) since they are non-zero at the boundary. Following 
H|, pH , it is convenient to project these functions on the space f/^J (f2; K 3 ), by subtracting the harmonic 
function on £1 with the same boundary data. Let ip : fl — > M 3 be the unique solution of the problem 



(21) 



= in n } 
cp = 8 on dfl, 



and set P8 = 8 — ip. We will often omit the dependence of ip on the parameters a, A, R, as for 8. 
From (|18|) and some standard computations it is easy to find that, in the case R = Id 



(22) 
and 



V = (jMZ, a) + o(A- x ), ~/i 2 (£, a) + oiX- 1 ), 1 - J^ 3 (£, a) + o(\- 2 ^j 



£^-l hl{ £,a)+o(\-i) 

(23) (*-¥>)=| ^^-ih 2 ^a)+o{\^) 

i+TFop +^fe3(Ca)+o(A- 2 ) 

where /ii and are defined in (|l^), and where /13 is the solution of 

(24) (A 6 h 3 (a,Z)=0 in ft, 

The quantities o(A _1 ) and o(A -2 ) in formulas (|22| ) and ( p3| ) denote functions which C fc (ft)-norm, for any 
k 6 N, is of order o(A~ 1 ) and o(A~ 2 ) respectively. 

We collect some further estimates, whose proof are trivial, and which we will use later. Given a fixed 
positive constant r < for A sufficiently large there holds 

(25) jf < f ; jf |W| < C^; J b \VP8\ < C^f; jf |,|» < C^; 

(26) + |V*|(sb) < T , Vxef7\B r ; |<p - (0, 0, l)|(a) + \Vip\(x) < T , Vx e 0. 

A A 

In the following, for brevity of notation, the constant C will be allowed to vary from formula to formula 
and from line to line. 

For k > 1 and for i,j € {1, . . . , fc}, i 7^ j, we will use the following notation 

e(e, A 1( . . . A fc ) = 0{e 2 ) + ]T 0(e\^\ log A,|) + ]T O (-L 



1,3 

e(e, A,) = o(e 2 ) + o^A' 1 ); e(A,, A,) = O ^(log A 4 + log A,-) ^ + ^ + + 



e( £ ,A 1 ,...,A fe )= ( e 2 ) + ^ ( £ A- 1 ) + £ e(A*,Ai) + ]T O (--^— 



We will often make use of the identity 

(1 + ICI 2 ) 3 

which is immediate to verify (the integrand is indeed the third component of A8 up to a constant). 



(27) / 7^^3=0, 
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3 A non-existence result via a Pohozahev-type identity 

In this section we prove a Pohozaev-type identity for the ff-surface equation. The proof is elementary 
and extends a previous result of Wente, see | 5jJ . 

Proposition 3.1 Let fl CM 2 be a smooth bounded and simply- connected domain, and let v G C 2 (f7; 
be a solution of 



(28) 



Av = H(£, v, Wv)v x A v y , in 57, 
v = on dil, 



for some continuous function if (£, Vv). Then v = in Q. 

Proof. We assume first that the domain 17 is the unit disk D. In the spirit of the Pohozahev identity, 
we consider the quantity X}j=i(£ ' Vi>j)Aui and integrate on D. We claim that 



(29) X^-VviXv* A «„),=(), Z=(x,y). 

i=l 

Once ( p9| ) is proved, we have 

3 

(30) ■ V ^)A^ = 0. 



i=l 



Integrating ( |30| ) over D and taking into account that the dimension is 2, we find 



3v ) I.^JdD \d> 

where v denoted the exterior unit normal to 3D. As a consequence we have |^ = on 3D. Thus, 
extending v to zero on the complement of D and also extending H continuously outside D we obtain a 
C 1 solution of Au = i?(£, k, Vu)?; x A v v in R 2 . Hence, applying Theorem 1 in Jl6| we obtain v = in D. 
Let us now verify (^9|): using simple computation we find 

3 

y^(£ • Vui)(^ A V y )t = X [{Vl) x (v2)x(v3,)y ~ (Vl )x {Vz) x (v 2 )y + {V2)x{vi) y (v 3 ,) x 
i=l 

- {V2)x(vi) x (v3)y + MxMxMy - (v 3 ) x (v 2 ) x (v 1 ) y ] 

+ V [{Vl) x {v2)x{vz)y - {Vl)x{vz)x{v2)y + (v 2 )x(vi)y(v3)x 

~ (V2)x(vi)x(v 3 )y + (V3)x(Vl)x(V2) y ~ (v 3 )x(v2)x(vi)y] = 0. 

This concludes the proof in the case fi = D. For the general case of a simply-connected domain, it is 
sufficient to use the Riemann Mapping Theorem and the transformation rule of (|2^) under conformal 
mappings. We recall that for 17 smooth, the Riemann map is also smooth up to the boundary, see Q. ■ 



4 The finite-dimensional reduction 

In this section we show how problem (^) can be reduced to a finite-dimensional one for small values of 
e. The starting point is the following Proposition, proven in [^lf (Lemma 5.5) using an isoperimetric 
inequality. We give an alternative proof in the Appendix, using the stereographic projection and shifting 
the problem from K 2 to S 2 . We believe that our proof could be naturally extended to the case of higher 
degree bubbles. 
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Proposition 4.1 There exists a constant Cq > such that 

l"(R5 a ,x)[R5 a ,x,R6 a ,x] < -C ||P<5 a , A || 2 , for all R6 a , x e Z. 

and 

T(R5 a , x )[v,v] > C \\\7v\\ 2 L2{R2) , for all RS a , x eZ and v ± (T RSa X Z ® {tR8 aX }t) 
In particular, the equation I (RS a x)[v] — implies v — c € Ti?<5 a X Z for some c€K 3 . 



We are going to consider now problem (15) on the domain f2. Given C > we set 
(31) 

Z = \ J2 PS PiAi : dist(pi, dQ) > C~\ dist(pi,pj) > C' 1 Vz 7^ j, X t e e [C~\ C], R, e 50(3) 



Proposition 4T asserts that the manifold Z, see (|l6| ) is non-degenerate for the functional I module 
translations. As a consequence, it is easy to extremize I e in the direction perpendicular to Z. This is 
stated in the following Proposition 4.3, in the same spirit as Q. We need first a preliminary Lemma (see 
also Q], Proposition 3.1). 

Lemma 4.2 Let k G N, C > 0, and let Z be as in (|3l|). Then there exists a positive constant C such 
that 

if v€H^(Q,),v±T z Z,v±P8iVi, then I"(z)[v, v] > C _1 ||t;|||-i {n) . 



Proof. For i = 1, ...,k, let Pj be the ball of radius around p, L , and let also Bi be the ball of 
radius around pi. Let us denote by Pi the orthogonal projection of Hq (£1) onto Hy{Bi), and for any 
v e Hq (fi) set «i = u — X)j=i It follows immediately that 



(32) 



ll w llH Mn)=Ell P ^ll^(n) + ll^llH Mn)- 



From standard regularity results, since the function v\ is harmonic in each Bi, and since it coincides with 
v on each dBi, there holds 

(33) IKIIc^s,) ^ C'lkllffoHo)' for alH = l,...,/c, 

where C is a constant independent of v. Since v is orthogonal to PSi, from (Eq) we deduce 



(P^, P6i) = \v-vi 



£i>,p<u = -(^^)+E°(f) IMk(n)- 



To evaluate the scalar product (vi,PSi) = J n Vui • VP8i, we divide the integral in the regions Bi and 
il\Bi. We have clearly Jo^s. V«i - VP5i = 0(X^ 1 )\\v\\ H i^-iy On the other hand, using (|25| ) and (|33|) wc 
find 

'log A, 



Vui • VP<5, 



Using these formulas and ( p6| ) we obtain 

'logA-i 



(34) 



A, 



+E° ( T I INI^SC") < Ce|log£| IM|j*j(n)- 
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In the same way as (p4|), using the explicit expression of the function Si and taking the scalar product 
of v with J^-, and J^, one finds ||IIjPji;|| < Ce\ log£|||i;|| H i( n ), where H denotes the orthogonal 

projection onto the space spanned by Si, J^f- and The functional I' E '(z) is given by 



I"(z)[v, v] = / Vv ■ Vv - 2 / z ■ (v x A v y + v x A v y ) + 2e / v ■ (g x Av y + v x A g y ) , v,v e H^fl). 
Jn Jn Jn 



It follows easily from the expression of /" and from Proposition 4.1 that 

(35) I' e \z)[PiV,Piv] > Co\\Pivf Hm - C £ |loge| \\v\\ B i (a) . 

For an arbitrary function v there holds 

k k 

(36) I'J(z)[v,v] = J2K\z)[PiV,Piv]+I'J(z)[vi,v 1 ] + 2Y,K\z)[PivM 

i=l i=l 
From the orthogonality of PiV and it follows that 

I'J(z)[PiV, Vl ] = -2 / z-((P l v) x A(v 1 ) y + (v 1 ) x A(P l v)y) + 2e P l v(g x A(v 1 ) y + (v 1 ) x Ag y ) 
Jn Jn 

= o(^J IzUVPiwIIVvxl) +0{e)\\v\\ 2 Hh{ny 
Dividing again the integral into the regions Bi and Q\Bi we deduce 

(37) rjizWv,*] = o IHI| i(n) + E° (x:) Mkw + °( £ )NIW 

Similarly, we obtain 

(38) i'M[viM = IKII^ (n) + o (^) Nl^ (n) + E° (xt) + 0(^)11^11^^)- 



From fl32|), (g5j), fl36|), (gTj) and (|38D the Lemma follows 



Proposition 4.3 Let C be a fixed positive constant, let k £ N, let e > 0, and Z be defined as above. 
Then, if e is sufficiently small, for every z £ Z there exist a function w £ (z) € ff 1 (f2; M 3 ) and C > wii/i 
£/ie following properties 

i) w e (z) is orthogonal to T Z Z, for all z € Z; 

Jg(^ + w e (z)) £ T Z Z , for all z £ Z; 

in) \\w e {z)\\ < C\\I' e {z)\\, for all z £ Z; 

By i) and ii), the manifold 

Z £ = {z + w £ (z) : z £ Z} 
is a natural constraint for I' e . Namely if u £ Z e and I' E \z E (u) = 0, then I' e (u) = 0. 



Proof. Given Proposition 4.1, the arguments are quite standard. For convenience, we give a brief 
sketch in the case k = 1. In the proof, we simply write 8 for R5 a \. 



Let us define F £ : Z x H^^R 3 ) x T Z Z -> ff^fijR 3 ) x K by setting 

F £ (z,w,q) -- 



I' e (z + w) — q 
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With this notation, the unknown (w, q) — (w e , V e {z + w E (z))) can be implicitly defined via the equation 
F E (z,w,q) = (0,0). Setting G E (z,w,q) = F E (z,w,q) - d( Wtq )F E (z,0,0)[(w, q)] we have that 

F £ {z,w iq )=0 d (Wiq) F e {z,Q,0)[(w,q)] + G e (z,w,q) = O. 

Reasoning as in Q, using Lemma 4.2 one can prove that d( Wiq )F E (z, 0, 0) is uniformly invertible for z G Z 
and e sufficiently small. Hence we can write 

F E (z, w,q) = & (w, q) = W E {z, w, q) := - (d (Wiq) F E {z, 0, O))" 1 [F E {z, 0, 0) + Q E (z, w, q)] , 
where 

Qe(z,w,q) = F E (z,w,q) — F E (z, 0, 0) - d {w>q) F E (z, 0, 0) [(w, q)] . 
It is also standard to prove that Q E (z, w, q) satisfies 



(39) 



\\Q E (z,w,q)\\ < CIK^,^)^ 

\\Qe(z,w,q) - Q E (z,w,q)\\ < C (\\(w, q)\\ + \\(w,q)\\) \\(w,q) - {w,q)\\ 



where \\(w, q)\\ and where C — C(£l,g,C) is a constant depending on fi, g, C, and independent of 
z G Z and e. Using ( |39| ) it is possible to prove that the function W E is a contraction in a ball of radius 
C\\F E {z, 0, 0)|| for some positive constant (7(0, g, C). Since ||F e (z, 0, 0)|| < C||i^(z)|| for some constant C, 
the conclusion follows. ■ 

We estimate now the quantity || PSi)\\ in order to control the norm of w E (z), see Hi) in Proposition 



Lemma 4.4 Let C be a fixed positive constant, let k £ N, let e > 0, and let Z be as in @. Then th 
holds 

||-f e ( z )ll — &( e i ^1) • • • i ^k)i for e sufficiently small and for all z G Z, 
where e(e, Ai, . . . , Xk) is defined in Section^. 

Proof. Let v S -ffQ(£l;R 3 ) and z S Z. Using integration by parts we deduce easily 

(40) l' E {z)[v] = / Vz ■ Vw + 2 / v ■ {z x A z y ) + 2e / z ■ (v x A g y + g x A v y ) + 2e 2 v (g x A &,). 

From Holder's inequality we get 

z(v x Ag y + g x A v y ) < 



ere 



HffllcHO) / N 2 HI, forallveiJ 1 ^ 3 ). 



From (p5|), ( p6| ) it is easy to check that (J Q |z| 2 ) 2 < CJ2i=i x"l log 1 2 , hence we have 



e I z ■ (v x A g y + g x A v y ) 
n 



< f |logAi|*||«||, for all v e H^(Sl; 



On the other hand, it is also immediate to verify the inequality 



< Ce 2 ||«||, for all v e H^(n) and all z e Z. 



It remains to estimate the first two terms in (f40[). Writing for brevity Si — Ri5p u \^ we have also 



Vz ■ Vu + 2 / v ■ (z x A z y ) 



J a v ■ [-a (x; ps.) + 2 ( (£ p* ) x a (x: pa,) t 
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Using the equation ASi — 2((5i) x A (Si) y ) and the fact that ASi — A(PSi), the above quantity becomes 

2 J a v - (J2 P6i ) x /K (X> P5i ) y ~^ 6i)x/K ^ Si)i 



which can be written as 



(41) 2 V / v ■ ((PSi) x A (PSj)y) I v ■ [{PSi - 6i) x A (S t ) y + (Si) x A (PS, - 

i7 tj J n i Jo 



Let us estimate first the term J n v ■ (PSi) x A (P5j) y . Let 7 < |C 1 (recall the definition of Z) be a fixed 
positive number and divide the integral in the three regions P 7 (pi), B 7 (pj) and O \ (B 1 (pi) U B 7 (pj)). 
Integrating by parts on the balls B 1 (p i ) and B 7 (pj), the quantity J„ u • (PSi) x A (PSj) y becomes 



Us • (PS t A (P^) B ) / v (PS, A (P*,-)x») + / v- (PSi A (Po^K 

B T ( Pi ) JB y ( Pi ) JdB^(p z ) 



v x ■ {PSi A (PSj)y) - / « • (PSi A (P^-)^) + / « • (PSi A (P^K 

v((PS i ) x A(PS j ) y ). 



'a\(B 7 ( Pi )uB 7 fe)) 

Hence, since Si and its derivatives are of order A^ 1 in f2 \ B 7 (pi) one finds 



V ■ ((P5i) X A (PS^y] 



< C 



1 



Using similar estimates we find that the whole expression in ( ]4l"| ) is of order e 2 . So we obtain the 
conclusion. ■ 



Lemma 4.5 Let C > and /ei Z oe as m (|3l|). TTien i/iere aoZrfs 

A- 



(9- 



9A, 



/or every i = 1, . . . ,k. 

Proof. From ( f4(i| ) and some integration by parts it follows that 

l"(z)[v, v}= / Vv ■ Vw + 2 / u • (z x A % + A z y ) + 2s / v • (g x Av y +v x A g y ) , 
Jn Jn Jn 

where v,v are arbitrary functions in iJQ(J7,R 3 ). We choose now v = J^-, and we let v be an arbitrary 
test function. We have clearly = *k _ where <ft is the function in ( |2l| ) corresponding to <5,;. 
From the estimates in (|26|) and from the explicit expression of jP 1 we find 



(42) 



dz 

-q^ ■ (g x A v y + v x A 5 y ) 



<C-|H|+Ce 



dpi 



\v\\ < Ce\\v\\. 



L 2 (n) 



Turning to the remaining two terms, we have 



„ „ dz „ f dz . , /" _ _ 9& „ /" cW, , . . . , . . . 

Vv • V— + 2 / — — (z x Av y +v x A z y ) = / Vv ■ V— + 2 —■ ((5 l ) x A v y + v x A (Si) y ) 
n dpi J n d Pl J n dp, J n d Pl 



(43) +2J2 7 r 1 -((Sj)^v y + v x A(S J ) y )+0(\^)\\v\\ 
&i J * dp > 
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Integrating by parts and using the fact that is of order A, ; 1 on dfl we find 



(44) 



Jn d Pi J n d Pl J 9n dv \d Pi 



0(X- x )\\v\\ 



since ^ satisfies (|o|). To estimate J n • ((^j^ A u a + v x A (Sj) y ) for j ^ i, we proceed as follows. Let 
Bi and Bj denote the balls of radius centered at pi and pj respectively. By the definition of Z these 
two balls are disjoint and moreover, by (^6|) V5i and V5j are of order A^ 1 and Xj 1 respectively outside 



Bj and B, . Hence we have 



(45) 



n d Pi 



{(Sj) x Av y +v x A(6j) y ) 



< C 



1 1 



Xi Xj XjXj 



IN 



Hence (42)-(45) imply I'J(z) -j^- < C \ e + J2j=i Xj 1 j ■ The remaining part of the statement follows 
from similar arguments. ■ 



From Proposition 4.3, critical points of I £ restricted to Z £ are true critical points of I £ . We define 
I £ : Z — > K as I £ (z) — I £ (z + w e (z)). We now analyze the reduced functional I e . 



Proposition 4.6 Let C > 0, let Z be as in ( |31| ) and let w e (z) be as in Proposition (.o. Then we have 



(46) 



I e {z)-I E {z) <C~e 2 (e,\ l ,...,Xk), Vz e Z; 



Moreover, for all z = RiP Pi ^x i € Z there holds 

(47) 



dhjz) _ dl e {z) 
dpi dpi 
_ dl c (z) 



dIAz) _ dh{z) 
dRi dRi 



< [e + j: j Xj 1 )e(e,X 1 ,...,X k y ) 
<i i (e + j: j Xj 1 )e(e,X 1 ,...,X k ) 



Proof. We have 

I e (z) - I e (z) = I s (z + w) - I e (z) = f I' e (z + sw)[w]ds = I' e (z)[w]+ [ (I' E (z + sw)- I' e (z))[w]ds. 

Jo Jo 

Since the functional I'J is locally bounded, we have the following estimate 



\I' e {z + sw)-I' e {z)\ < 



I"(z + tsw)[w]dt 



< sup \\I" (z + tsw)\\\\w\\ < C\\w\\ 

t,se[o,i] 



for some fixed con stan t C depending on f2, ||g||c 2 (<9£i) and the above constant C. Using the last three 
equations, Lemma 4.4 and the property Hi) in Proposition 4.3 we find 

le(z) - I s {z) < \\I' s (z)\\ \\w £ (z)\\ + C\\w £ (z)\\ 2 < i 2 (e, Ax, ... , A fe ). 



This concludes the proof of (46). We just sketch the proof of (47). Differentiating the equation 
F £ (z, w, q) — with respect to pi we obtain 

dF £ dz dF £ d(w,q) „ dz dF £ d(w,q) 

dz dpi d(w,q) dpi dpi d(w,q) dp t 

Similarly as before, one finds that d^ w ^F £ is uniformly invertible, and hence 



(48) 



dw 


< c 


dpi 





< c 



T»( \ 9Z 
lAZ) dp~ 



C\\w £ (z)\\ 



d?: 



dpt 
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where we have used the fact that I'J is locally Lipschitz. We have 

dl e (z) dl £ (z) ra ~ in 



d Pi d Pi : lAz) 
(49) + I'M 

Equation (^9|) implies 

di e (z) di e (z) 



dz 

-K-,w e (z) 



+ I (I'J(z + sw e (z))-I'J(z)) 





dz 

■^—,w e (z) 

OPi 



dw 



dpi 



dpi dpi 



C 



dw 
dpi 



dz 



,w e {z) 



+ / (l' E '(z + sw E (z))-l'J(z)) 
Jo 



dw 
dpi 



w e (z) 



+ c\\i'M 



dpi_ 
dw 



dpi 



\w e {z)\\+C\\w e (z)\\ 2 
dw 





dz 




dw 




( 


dpi 


+ 


dpi 


) 



■C\\w e (z) 



dpi 



Then the estimate of ^§^- - ^§^- in @ follows from Lemma fH], (||) and Lemma fO|. The remaining 
part of (MTu follows from similar estimates. ■ 



5 The expansion for one bubble 

In this section we compute the expansion of I e (z), with z G Z, for e small and in the case k = 1. This 
is essentially performed in p0[-|2l|, in order to construct blowing-up solutions of (^), and in order to 
characterize the mountain-pass solutions in the limit s — > 0. We derive the expansion here, in a form which 
is useful for us in the expansion for multiple bubbles in Section Let us first introduce some notation. 
We recall that g : ft — > R 3 denotes the solution of (||), and letting R e SO(3), define dug : ft — > K by 

d d 
dRg(a) = —{Rog) 1 (a) + — (Ro g) 2 (a), x e CI. 

For a fixed boundary datum g, we are interested in expanding the functional I e (PS) as a function of the 
parameters a, A, R and s. We have the following Proposition. 

Proposition 5.1 Let C > be fixed, and let a, A, R be such that PR5 a ,x G Z . Then, setting 

Id 2 



An 



(1 + KI 2 ) 3 



r it „\:.o.\.Jt} = s.S„ ( —H( a ) - jd R -ig(a) 



there holds 



dI £ (PR6a,. 



I E (PR5 a , x ) = -A + F n Ja, A, R) + o(e 2 ) + e(e, A); 
y 



dF, 



+e(e, A); 



dI £ {PR5 a , x ) _ dF n , g , e(e,A) dI £ (PR8 a . x ) _ dF n<g 



d\ 



dX 



x 



dR 



dR 



+e(e, A), 



da da 
where e(e, A) is defined in Section |j. 

Proof. We assume that R — Id, and we write 8 for 8 a \. Let also ip be the solution of (|2l|). We have 
I e (PS) = ~ / |V<5| 2 + i f |V^| 2 - / V^-W+^ f(6-<p)-((6-tp) x A(8- 



(50) 



e I (8 -tp) -((8- tp) x Ag y +g x A(S- tp) y ) + 2e / (8 - (p) ■ {g x A g y ). 
in Jn 



Integrating by parts we can write 

(51) \ f |V<5| 2 + 1 f |V^| 2 - f Vp-V8= f (cp — 5) ■ (8 X Ad y ). 

* Jn z Jn Jn Jn 
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We expand first the expression in (pl| ). Let us evaluate the z-component in the scalar product of the 
integral on the right-hand side in (|5l[). From formulas (|l8|) and (23) we deduce 



(<p-S) 3 (5 x AS y ) 3 



Integrating on f2 we get 



, 2 h i ^ A 2 jC-q| 2 - 1 
l-^ 3 (a,a)- A2|e _ a|2 + i 



+ o{\- 2 ) 4A 



■ l-A 2 |£-a| 2 
(A 2 |£-a| 2 + l) 3 



-8X 2 



A 2 |£-a| 2 + l A 2 



— h 3 {a,a) +o(A 2 ) « 



(52) / ( v ,-«5) 3 (4a^) 3 = -8A 2 



2 , l-A 2 |e-a| 2 
^ /n(A 2 k-a| 2 + l) 4 



1 



Using a change of variable we obtain 

'/ 

Jn 

(53) 

and also 



A 2 |£-a| 2 + l A 2 

l-A 2 |£-a| 2 



- rnh 3 (a,a) + o(A 3 ) 



= A 



, (A 2 |e-a| 2 + l)i 



R 2 \0 



( A 2|£_ a |2 + 1) 3- 

l-A 2 |g-a| 2 
(A 2 |e-a| 2 + l)3- 

1-A 2 |C-«| 2 
(A 2 |£-a| 2 + l) 4 



I + 0(A- 4 ), 



l-A 2 |e-a| 2 
(A 2 |£-a| 2 + l) 3 



(IC| 2 + 1) 

l-A 2 |£-a| 2 



(54) 



(A 2 |£-a| 2 + l) 3 
1 /• 1-ICI 2 



1- A 2 |£-a| 



(A 2 |£-a| 2 + l) 3 
+ 0(A- 4 ) = 0(A- 4 ). 



a^ j R2 dei 2 + 1)3 

The last identity follows from (p7|). In conclusion, from (|52|), ( |53"1 ) and (|54| ) we get 

IC| 2 -1 



(55) 



(<£ - 5)3(4 A Sy) 3 



* (l^l 2 + I) 4 



o(\- 2 ). 



We consider now the x and y components of the integral on the right-hand side in (|5l]). We have, using 
© and © 



We can write 
A 4 



(99 - S)i(S x A S y )i = 
(x - ai) 2 



2, x 2A(a;-ai) 
A fel( ^ a) " l + A 2 |^-a| 2 +0(A } 



(l + A 2 |£-a| 2 )^ 



= A 4 



(x - ai) 2 



8A 3 (a; - ai) 
(l + A 2 |£-a| 2 ) 3 ' 

(.t - ai) 2 



(l + A 2 |C-a| 2 )4 Mo (l + A2|e-a| 2 r 

<- 2 

+ 0(A- 4 ). 



(1 + iei 2 ) 4 



From the smoothness of hi we have also 

(hi(£, a) - hi(a, a) — — a) ■ Vhi(a, a)) 

As a consequence we deduce 
(x - ai) 



x — a± 



hi(x, a) 



(l + A 2 |C-a| 2 ) 3 



h\(a, a) 



(l + A 2 |C-a| 2 ) 3 
(x - ai) 



< O 



K-4 



(l + A 2 |£-a| 2 ) 3 



(l + A 2 |^-a| 2 )3 



+ /((£- a)- V/ii (<*>»)) 



(x - ai) 



(l + A 2 |£-a| 2 ) 3 

d x hi(a,a) 



(1 + A 2 |e — a| 2 ) 3 / A 4 



(1 + ICI 2 ) 3 



+ o(A- 4 ). 
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From the last equation we deduce 

(56) Jjfp - 5U5 X A Sy), = 16 £ (1+ ^ |2)4 - 16^ d xhl (a, a) £ (1 + ^ |2)3 + o^ 2 ). 
In the same way we obtain 

(57) JW - S U S, A « 2 = 16 / ri ?T ^1 FF - l 6 i W „, .) £ IT ^ F + o(A-). 
From (||), © and (||) it follows that 



(<p-6)-(5 x A5 y ) = 8/ 7 J^—L +32 



It is standard to check that 



(i + iei 2 ) 4 u (i + m 2 ) 4 

'4(1mp?)(I^K")-<»-')' 



(i + iei 2 ) 4 7 R2 (i + ICI 2 ) 4 3 

Hence from the last two equations we find 

(58) j^: - 5) ■ (S x A 5 y ) = - J^ + ^) (,. a ) • >! A ! ». 



We turn now to the fourth term in (p0|). We have clearly 

(5 - (p) ■ {{6 - y>) x A (S - if) y ) = (5 - ip) ■ (S x A S v ) - / (<5 - tp) ■ (S x A ip y + <p x A <5 a ). 
2 Jo Jo 

Let us consider the term 6 X A ip y . Using the above formulae we deduce 

(X a,„\ 8A(x-ai)(y-a 2 ) ( dh 3 nn -iC\ 8\{x_ai) ( dh 2 , nn -i,\ . 

8A(a?-ai) /9/n u \ 4 1 + A 2 (y 2 - X 2 ) /3/i 3 

(59) (U ^ = FCTl¥ +0(A ) J + A (l + A 2 |"-aP)H ^ + ° (A V ; 
. , . l + \ 2 ((y-a 2 ) 2 -(x-ai) 2 ) (0h 2 \ SX^x - ai )(y ~ a 2 ) (dhx x A 

(4A ^ )3 = 4 (i + A 2 ie-ap) 2 l^ + 0(A V + (i + A 2 ie-a| 2 ) 2 ^ + ° (A V' 

and also 

8A(y-q 2 ) /dfo 2 \ 4 1 + A 2 ((x - a,f (y a 2 f) ( dh 3 
( ^ A ^ )l= (l + A 2 |C-a| 2 ) 2 +0(A V + A (l + A 2 |C-a| 2 ) 2 [l* + ° {X ] 

(an\ (,n am _ 8A(g-ai)(y-Q2) (dh 3 nfy -x^\ 8\(y - a 2 ) ( dhx nn -i, \ , 

(60) ((fx A 6 y )2 — (1 + A2|e _ a|2)2 ^ + 0(A )J- (1 + A 2| e _ a | 2)2 ^ + ))' 

, . r \ , l + A 2 (( 3 ;-a 1 ) 2 -( j/ -a 2 ) 2 ) (dhx , „ r ._ 1 A , 8A 2 (x - oQCg - a 2 ) /0fe 2 x 
C^ A ^)3 = 4 (l + A 2 |e-a| 2 ) 2 ^ + ° (A V (i + A 2 |C- fl | 2 ) 2 ^ + ° (A } 

The functions Ip 1 in the last formulas, as before, are evaluated at the point (a, a). Using (123), (pit) and 



0|) we find 



(01 ) / (8-<p). (S x A Vy ) ~ -32i p^jSjs) ^(«. ») + 
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and similarly 



(8-<p)- { Vx A S v ) ~ -32^ 



— (a,a)+o(X 2 ), 



Since A <p v )i\ , A <p y )i\ < 0(A -3 ) and A <^)3| < 0(A~ 2 ), one can check that 

(d-<p)- (cp x A <^ a ) = o(A~ 2 ). 



Let us now turn to the fifth term in (|50|). The quantity (# — 95) • ((5 — ip) x A g y ) can be estimated as 



2A(x-oi) 2 

— T"-l 



1 + A 2 |e — a| 2 A 



4A 3 (x-g 1 )(y- fl2 ) 2 , 
- (1 + A 2| C _ a | 2)2 - ^ ) (33), 



4A 2 (x- ai ) 2 , . 

+ — <9^3 ) (.92), 



(l + A 2 |£-a| 2 ) 2 A 2 

2A(y-a 2 ) _ 2 
1 + A 2 |e — a| 2 A 



4A 2 (x- ai ) 2 

+ — <9^3 ) 



(l + A 2 |£-a| 2 ) 2 A 2 



, l + A 2 ((y~a 2 ) 2 -(x- ai ) 2 ) 2 
2A (l + A 2 |£-a| 2 ) 2 A ' (53)a 



A 2 



1 + A 2 |e - a| 2 



2A 



l + A 2 ((y- a2 ) 2 -(x- ai ) 2 ) 



(l + A 2 |£-a| 2 ) 2 



+ ^ (i + A 2 |£-a| 2 ) 2 +a^V (51 ^ 
Integrating on 51 and reasoning as before we get 



(62) 

and similarly 
(63) 
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(6-<p)- ((6-<p) x A g y ) = — - 



(S-cp)- (g x A (S-tp)y) 



16 

T 



(i + ICI 2 ) 3 



a + iei 2 ) 3 



A 



d x hi (g 2 ) y 



(g 2 ) y + 0(\- 2 ), 



( 9l ) x + o(\- 2 ). 



Finally, the last term in ( |50| ) is easily seen to be of order o(e 2 ). This concludes the proof in the case of 
R = Id. For a generic rotation R it is sufficient, by invariance, to consider the boundary datum R~ 1 g 
and to substitute (gi) x + (52), with d R -ig. ■ 



Remark 5.2 Propositions ^.6 and 5.1 allow us to find critical points of I £ extremizing the reduced func- 
tional I £ on Z . Differentiating with respect to R, X, a we get 



(64) 



■^d R -ig = 0: 2H(a)-e\d R -ig = 0; ^V#(a) - jVd R -ig = 0. 



Using the second and third equations in (Q) we deduce 

V log H(a) = 2V log d R -ig(a). 

The extremization with respect to R (the first equation in (p4[) ) is performed in MJj Lemma 5.4 and |§ 
Lemma 3.1.2 and, requiring d R -ig to be positive (from the second equation in (|64|) ) yields 



(65) 



d_ 
dR 



d R -ig = 



d R -ig= (|V 5 | 2 ±2| 5;c Ag y \Y 
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Hence, under the conditions V<? ^ 0, the extremization in ( |64|) becomes 

ed R -ig(a) oR H 

In particular the mountain-pass solution of (||) (s ee %2l\] ) has minimal energy on Z , and one has to chose 
the + sign in (|6^)-(|66|), and the last function in (|66|) is maximized on St. 



Remark 5.3 (a) We point out that the expansions in ( pl| ) and (£>8|) yield 

\\PS\\ 2 HH u) = ¥\\l-^H(a) + o(X- 2 ). 

Since the norm of P5, being the projection of 5, is smaller than the norm of S, the above formula implies 
H > on SI. A little more calculation shows that indeed H > 0, as proved in jfefy. 

(b) The functions hi, hi, and /13 defined in (|l2|), ( p4| ) are related to the boundary values of 8. Since 
hi and hi are of order A -1 , while hi is of order A -2 , /13 appears in the expansion only as a lower order 
term. 



6 The role of the Robin function 

In this Section we investigate the relation between the Robin function H and the function H defined in 
(fl4|). Since H consists of second derivatives of the regular part of the Green's function, while the Robin 
function involves the regular part itself, we need to use global arguments, based on the Riemann mapping 
theorem. See || for some properties of the Robin function. 

6.1 Simply connected domains 



In this subsection we prove the first assertion of Theorem 1.1 



Proposition 6.1 Let SI C M 2 be a smooth simply connected domain. Then, if f : SI — > D is a Riemann 
map, there holds 

(67) e^)= 1« ■ H(a)=2 ^ 



(l-|/'(«)l 2 ) 2 ' (l-|/'(«)l 2 ) 2 ' 

PROOF. The first part of the statement is well-known, see e.g. jfjj, Table 2. Letting Gd(el,C) denote 
the Green's function for D, and letting <p : D — > M being any smooth function with compact support, we 
have 



<p(a')= f G D (a',0A<p(0d£, a 1 ED. 

J D 



ID 

Let %/} : SI — > M be defined by ip — f f ■ We have A^(£) = jjtj^i A</?(£) and hence, letting G(a,£) be 
the Green's function for SI and using a change of variables we get 

1>(a) = <p{f(a)) = <p(a') = f G D (a', 0A V (0« = / G D (f(a),f(0)^H, 

Jd Jn 

where a' = f(a). Hence, from the explicit expression of Gd it turns out that 

G(a, z) = G D (f(a),f(z)) = -\ log j f{z] ^ - ; a, z e SI, 

2 |i-/(2)/(a)| 2 



where we have identified £ with the point z in the complex plane. It follows that 

a. z G SI. 



H{a,z) = \ 



logl.-a^-log.^-^l 2 



f(z)f(aW 
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The last expression can be rewritten as 



H(a,z) 



log\l-f(z)f(a)\ 2 -lo i 



\m-f(a)\* 



\z - a| 5 



s,z6fl. 



In particular, taking the limit z — > a, we deduce immediately the first equality in (|6 
Using complex notation, we have 



_d__(d_ d_\ d _ ,( d d\ ^__(^_ &\. A- Y— - — 

dai \da da J da 2 \da da J ' 9a; \ dz dz J ' 9y \ 9z 9z 

It follows that 
(68) 



d d d d 3 3 „ d d m 9 9 
1 = 2 + 2— — = 4Re— — := L. 

dx da\ dy 9a 2 dz da dz da dz da 



To derive the expression of H , recall (|13[), we apply L to H(a, z) and evaluate at z = a. We have, still in 
complex notation 

8 = 1 7Mf(a) _ 19 /(,)-/(q) 

9a v ; 2i_/(2)/( fl ) 2 9a z-a 

When we apply the operator the second term vanishes and we get 



9 9 1 /WOO 

T^T^-n(a,z) = — — = • 

9z9a 2(i_/( 2 )/( a ))2 



Taking the real part we find 



9 9 

4Re— — i?(a,z) 
9z 9a 



(l-f(z)f(aW (l-f(a)f(zW 



Choosing z — a in the last formula, we obtain the second identity in (67). This concludes the proof. 



Remark 6.2 (a) The expression of H in ( p7\ ) does not depend on the choice of the conformal map f of 
Q into D. 

(b) From the explicit description of H in Proposition 6.1 we obtain H(a) — > +oo as a — > 957. This is 
true for any domain, as proved in $2l\j, Lemma 5. 7. 

(c) In the case of simply connected domains, the function H coincides with the square of the reciprocal 
of the conformal radius and the hyperbolic radius, see W, Definitions 1, 7 and Theorem 8. See also 



Remark 6.4- 

(e) Since every convex domain has a single conformal incenter, see JZ^ / Proposition 11, it follows 
that H possesses a unique critical point in this case. For a general simply connected domain H will have 
multiple critical points, see J7^ / page ^83. We also point out that, even if a conformal transformation of 
the domain affects the number of critical points of H , the topology of critical points at infinity (see MJ) 
at the first level of non- compactness should be an invariant. 



6.2 Multiply connected domains 

In this subsection we derive a general formula for H on multiply connected domains. This formula makes 
use of the covering map and the deck transformation. 

Let us recall that the Green's function in the unit disk with pole zq £ D is given by 



G D (z,z ) = - log 



1 - Z Z 



z, z £ D. 
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Let us pick a point w E fl and consider the Green's function for with pole at a. From ||, Theorem 4, 
one has 

Gn{a, w) = — log - k _ — , where f(z) = w, f{zk) = a for all k. 

k 

Thus the regular part Gci(a, w) is 



1 - z k z 



(69) H n (a, w) = log \f(z ) - f{z)\ - £ log \z k - z\ + ^log|l 



z fc z|. 



Now, as in the previous subsection, it is sufficient to apply the operator L defined in (p8|). The first two 
terms vanish when L is applied. To handle the third term, note that / _1 is a local diffeomorphism. So 
/ _1 is defined from / _1 : W — ► Vfc, where U is a neighborhood of a and Vfc is a neighborhood of Thus 



we have 

Similarly, there holds 
(70) 



dH _ dH dz k _ 1 dH 
da dz k da f'{z k ) dz k 



d 2 H 1 1 d 2 H 



dwda f'(z k ) fi( z ) dzdz k 
Hence, using (|6^) and ([70]) we find 

d 2 H , _ 1 1 



4Re 



2ReV L= 



A' 

Let Tfc be the Mobius (deck) transformation that maps Zo m t° z fc- We have 

f = f°T k f'(z ) = f(z k )U(z ). 

Using the last equation and factoring the term (1 — |zo| 2 ) 2 , we deduce 

[n) /H - e dwda lw = Wl ~ V(*>)l 2 (i-N 3 ) 2 ^ (i-z k z r ■ 

From ( f7l| ) we obtain immediately the following result. 

Proposition 6.3 Let il C R 2 fee a multiply connected domain, and let f : D — > il fee a conformal covering 
map. Given a £ tt, let {z k } k be the pre-image of the point a under the map f , and let T k : D —* D denote 
the deck transformation mapping zq into z k . Then there holds 



(79] fr (n) 1 v T k (z Q )(l~\z \ 2 ) 2 , T k (zo)(l-\z \ 2 ) 2 



i/'(zo)i 2 (i-koi 2 ) 2 ^ v (i-z k z f (i-^zo) 2 ; 

We note that when f^ 1 (a) = {z Q } we recover the formula for the simply connected domain. 
6.3 Some numerical computation 

In this subsection we prove that in general, for a multiply connected domain, the two functions H and 
2e 2H do not coincide (this is the case for simply connected domains, see Proposition 6.1). We consider 
in particular the case of an annulus of inner radius ^ and outer radius p, where p > 1. Our numerical 
computations show that t he c ritical points of these two functions do not coincide, hence we obtain the 
statement (b) in Theorem LI 
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For p > 1 we set 

1 



A p = l(x,y) : ^ < x 2 + y 2 < p 2 \ ■ S p = {(x,y) : -logp < x < logp} 



P 

It is clear that S p is a covering of A p through the exponential map. We also define a S C, h p : S p — > D 
and / p : D -> A p by 

a = -^-j ; Op u; ) = — — — r, ftp z =-log - , f = expoh p . 

2 log p e aw + 1 p a \ 1 — z J 

where w £ S p and z G D. Our aim is to compute formula J72] ) for this particular case. Fixing zo £ D, 
the points z^ and the corresponding points Wk = h~ 1 Zk are given by 

z fe = ? r ; w k = — log + 2fc7TZ. 

Using some elementary computations we obtain 

z + M k ( kir 2 

z k = — 77 — —', where Mk = tanh 



ZoM k + r \2\ogp 
If Tk denotes as before the deck transformation, then there holds 



zM fe + l' fcw (1 + M fc z) 2 



By symmetry, it is sufficient to compute ( |7l| ) for /(zo) real and positive. It is convenient to use the 
following parametrization for the points a € A p and zq € Z? 

a = logx; z = -itan(— - — logs], x € (- logp, logp). 

\41ogp J 

Using this notation, from equation ( |72|) we are left with 

1 W (1-M fc 2 )(l-N 2 ) 2 (l-Af fc 2 )(l-|*o| 2 ) 2 



^ Lf / (zn)l 2 (l-|z n | 2 ) 2 ^ 1(1 



|/'(z )| 2 (l - |z | 2 ) 2 V V (! + ^o) 2 (l - Zfezo) 2 (1 + M fc z ) 2 (l - z fe zo) 2 

From the above formuls it follows 

fr ( , 2 v (1 - Mj)(l - |zq| 2 ) 2 ((1 - |zq| 2 ) 2 - 4|z | 2 M 2 ) 

W |/'( Zo )| 2 (l-|z | 2 ) 2 ^ ((l-|zo| 2 ) 2 +4|z | 2 M 2 ) 2 

From [p| we have 



e 



2H{a,a) 



n 



|/'(z )| 2 (l-|z | 2 ) 2 



Zfe - Z 



1 - ZfeZ 



Using elementary computations it turns out that 



2 -i / 00 

(73) g(a) = R ,/ y2 7 \ 1 + 2^W(M) 

8(logp) 2 CO s 2 (^logx)N 2 V 



fc=i 



(74) 2e 2 "^) = l - , J] Z(k,x) 2 

8(logp) 2 C0S 2 f_?_ logx \ N 2 fe ^ 
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- ">" 4 (life 'oe 1 )) + 4 ta » tl (nfe l0 «*) 

In Figures 0^ we plot the functions 2e 2H and H (modulo the irrelevant factor 8 ^ gp yi ) for p = e and 

for p = e 3 5 . We note that, roughly, W(k,x) ~ e~ and Z(k,x) ~ 1 — e _T °ip so for small values of p 
the terms with k ^ are almost negligible. This accounts for the fact that for p — e the graphs are very 
similar, see Figure |l|, even on a fine scale, see Figure [l| For large values of p the difference between the 
two functions is mpre pronounced, see Figure ^. 

Remark 6.4 We recall that the harmonic and hyperbolic radii are defined by 

r har (0 = e- H ^; r hyp (f(z)) = \f'(z)\(l - |z| 2 ), 

where f : D — > Q denotes a conformal covering map. See ^ Definition 1 and page 15. Note that 
the function H, in the case of general non-simply connected domains, do not even coincide with rj~jL. 
However, in the case of small annuli, numerical computation show that rhar and rh yp are very close, see 
Figure 8. We also point out that the harmonic radius is related to the Bergman kernel, see J^/ Section 
8.1 



7 The expansion for multiple bubbles 

In this section we consider the case of multiple bubbles. We begin considering only two bubbles i?i<5 Q ,Ai 
and i?2^b,A 2 i which we denote for simplicity by 8\ and 62 respectively. We assume that R\ — Id, namely 
the first bubble is not rotated, and we simply write R for R 2 . 
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Figure 2: the functions 2e and H for p — e, detail 
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For C > 0, k = 2 and P8\ + P5 2 G our aim is to expand the functional I s (see (|Pj) on Z in terms 
of the parameters a — pi, b — p2, Ai, A2 and i?. In the following, for brevity, we set (see Section ||) 

(7 = 6- a, e(A x , A 2 ) = O ( (logAi +logA 2 ) ( i + i 



A^ A2 A^ A2 A1A2 
We recall the explicit form of the functional / e (w), for u G i?o (^5 

I e {u) = 



|Vu| + -/ u-(u x Au y ) + ej u ■ (u x A g y + g x A u y ) + 2e j u ■ (g x A g y ). 



7.1 Interaction with g 

We consider first the interaction term e J n u{u x A g y + g x A itj,) in (|^|), with u = P<ji + P5 2 . We recall 

that throughout this section we assume that C A^ 1 < e < CAr 1 for some fixed constant C and for 
2 = 1,2. We have 

i ' («x A 5y + fe A u y ) = [ {P5 1 + PS 2 ) ■ [((PSi) x + [P5 2 ) x ) Ag y + g x A {(PSi) y + [P5 2 ) y )\ 

Jn 



(75) = P5i- ((P5i) x A g y + g x A (P5i) y ) + P5 X - ((P5 2 ) x A g y + g x A (P&)») 

Jo Jo 

+ P8 2 - ((P6 2 ) x A fls -(- g x A (P* a )») + / • ({PSi) x A fls + & A (P*i)„). 
Jn Jn 

The first term in ( |75| ) has been estimated in Section [| formulas (|62|)-(p3|), and gives 



(76) 



e / P5i((P8 1 ) x Ag y +g x A{P5 1 ) y ) = -l§- 



Ai in (1 + ICI 2 ) 3 



d/dff(a) + O(A^). 



The third term in (|75j) can be estimated similarly, using the invariance of the problem under rotation, 
and gives 



(77) e f P5 2 ({P6 2 ) X A g y + g x A (P6 2 ) y ) = -16^ / 

Jo A 2 Jn 



o(l + |Cl 



2VS 



^- 1 . 9 (6) + 0(A 2 - 2 ). 



Let us compute now the remaining two terms in (|75|), starting from the fourth. We write a = {ai,a 2 ) 
and b = (b%, b 2 ). Up to an error of order e 2 , we have 



(78) (PS 1 ) X A. 



and 



4\ 3 1 (x-a 1 )(y-a 2 ) 2 F) h \ ( n ~\ ( 4Af(g-ai) 2 a , \ / \ \ 

' (i + A ^_ a |2)2 " Xl d * h 2J ^3)y ~ ^ (1+A 2| ? _ Q | 2) 2 + A^^/lsJ (<? 2 )y 



((ifx^-ah* + i? 9 *^) (ffl)y 



l+A^((y-a 2 )"-(x-a 1 )") 2 



P5, 



'ni ?"i2 r i3 

7"21 ^22 ^"23 
v?"31 ^32 ?"33 



2A 2 (;/-& 2 ) 2 L f\ 

l + A||C-f,p _ X7«2(A<U 

V Aj' 13 ~ T+A^-&F 



where are the entries of the matrix R. We are going to prove that 



(79) 



P5 2 {(PSi) x A g y + g x A (Pfc) v ) = O (e 2 \ logeQ . 
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If one uses ( f78| ) and (|7^), the integrals involved in the above expression are of the form 

(x - ai){y - a 2 )(x - &i) 3 / (x - ai)(y - a 2 )(y - b 2 ) 



(80) A?A 2 

(81) A? 



Af A 2 



q (1 + A?|C - a|2)^ (1 + AllC - &| 2 ) ^ (l + A?|e-a| 2 )'(l + A2|e-6P) 

- ai)(y - a 2 ) A? /" (x - ax)(y - a 2 ) Af f (x-a 1 )(y-a 2 ) 



n (1 + Af |C - a| 2 ) 2 (1 + A|1C - &| 2 ) ' A 2 (1 + Af|^ - a |2) 2 ' A 2 7 n (1 + A?|£ - a| 2 ) 



(82) 
(83)A? 



A 2 A 2 



(a; - ai)(a; - 61) 



A 2 A 2 



(x - at)(y - 6 2 ) 



a (1 + A 2 |C - a\'Y (1 + A||«e - 6| 2 ) ./n (l + Af |£ - a| 2 ) z (1 + A 2 |£ - 6| 2 ) ^ 

(x — A 2 /" (x — ai) Af /" (x — 



(l + A 2 |e-a| 2 ) 2 (l + A 2 |e-6| 2 )' A 2 7o (l + A 2 |C-a| 2 ) 2 ' A 2 Jq (1 + A 2 |£ — a| 2 ) 



(S4)A|A) ; l + A 2 ((y-a 2 ) 2 -(x-a 1 ) 2 )( a :-6 1 ) ; ^ f 1 + A?((y - a 2 ) 2 - (a; - ai ) 2 )(v/ ■ /, > ) 



(85) 



n (l + A 2 |e-a| 2 r(l + A 2 |C-6| 2 ) ./n (1 + A 2 |£ - a| 2 ) z (1 + A 2 |£ - &| 2 ) 

1 + \l((y - a 2 f -(x- a.f) Ai /" 1 + A 2 ((y - a 2 ) 2 - (s - ai) 2 ) . 



(86) ^ 



(87) 



Ai 



Ai /■ l_+A 2 (( 2/ -a 2 ) 2 -(x-a 1 ) 2 ) A 2 

2 "'SI 



(l + A 2 |^-a| 2 r 

1 r 1 



x — 61 



A, 



y - h 





(1 + 


A?ie- 


a| 2 ) 2 


A 2 




y-b 2 




Ai 


X (14 


All?- 


6| 2 )' 


A 2 




x — 61 




A? 


1(14 


All« - 


6| 2 )' 




/ 1 




1 




'n A 2 A 2 


; / 

Jn 


A? A 2 , 



A 2 7^(l + A 2 |?-6| 2 )' A 2 7 (l + A 2 |?~6| 2 )' 7n A 2 A 

The errors in the expressions of {P8i) x A 5 y and PS 2 are negligible with respect to the quantities listed 
in (|80|)-(|88|), hence it is sufficient to consider the above expressions. 

Estimate of (|80|). Using the rescaling £ 1— ► Ai(£ — a) and setting Oax,o = Ai (SI — a), we get 

{x - a\){y - a 2 ){x - b x ) _ A 2 



A?A 2 



a (1 + Af |£ - a| 2 ) 2 (1 + A 2 |e - 6| 2 ) A? ^ a (1 + |^ 2) 2 (l + - Al a| 2 ) 



x 2 y 



(89) 



A 2 

-7-O-X 



xy 



A i V.. (i + l£l 2 ) 2 (i + ^-AH 



Let us consider the first integral. We divide Q Xl ,a into the regions |£| < ^ an d |£| > if |£| < Ai2l 5 
then I? — Aiffj 2 > , so we have 



A2 
A? 



x 2 y 



i«i<^ (i + iei 2 ) 2 (i + f ic-w 



< 



CA 2 



CAi A2 

A 2 A 2 |a| 2 4,<^ (1 + |£| 2 ) 2 " AfAfFF 



< 



where C is a positive constant depending only on f2. When |£| > Al j g l then for Ai large there holds 

l$l 3 
(I+IfP? 



— < Ai 2 er | . As a consequence, using also the change of variables ^-(6 — Aicr) 1— > £, we deduce 



A2 
A? 



f a^y 

'\i\>^ (l + |C| 2 ) 2 (l + f ie-Aia|- 



< 



CA2 / Ai 



1 <c ,_logA 2 



A?k| VA 2 ; V , (l + |e| 2 ) " A 1 A 2 |a| 
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Turning to the second integral in the r.h.s. of (|S9|), we again divide the domain into two regions |£| < 



and |£| > ^f-. Reasoning as before we find 



A 2 



xy 



< C- 



A 2 (Tl 



ICI 2 < C K| logAx 



AiA||cr| 2 7^|<Aifi (l + |e| 2 ) 2 - M 2 AiA 2 



A 2 

AT ai 



Cl> 



< C 



A 2 <7i / A 



1 < ^l^il A 2 logA 2 



*2.<> 



(i + iei 2 )- H 2 a? 



Since in the definition of Z we assume dist(pi,pj) > C , \a\ is uniformly bounded from below, the last 
formulas imply 

(x - ai)(y - a 2 )(x - h) 



AfA 2 



la (1 + Af |e - «| 2 )" (1 + AHe - &| 2 ) 
The second expression in ( |so|) can be estimated in the same way. 



< C^llog el. 



Estimate of (pl|)-(|88|). We just treat some particular cases, since many terms are similar to each- 
other. First of all, all the terms for which the quantity (l + A 2 |£ — a| 2 ) 2 (l + A?,|£ — 6| 2 ) appears in the 
denominator can be treated as before. 

Next, we consider for example the second term in ( |S5| ) and the last term in (|6|). Using the changes 
of variable Ai(£ — a) i— > £ and A 2 (£ — b) i— > £ we find 



Ai 

A 2 Ja 



1 



l + Xl{{y-a 2 Y -(x-arf) 

(i + A 2 ie-«i 2 ) 2 

A 2 f y -b 2 



1 + y2 -j <C^<Cs^o g s\; 



AiA 2 y nXiiQ (1 + |£| 2 ) 2 " AiA 



Ai7o(l + Aiie-6| 2 ) A!A 2 Jn 



< C- 



1 + |C| 2 " AiA 2 



< Ce 2 . 



Conclusion. Using the estimates of ( |80[ ) and those of (|8l|)-(|88j), we obtain (|79|). In the same way, one 
can prove that 



(90) 



PSi((PS 2 ) x Ag v +g x A (PS 2 ) y ) = O (e 2 \ logeQ 



Hence, from equations (|76|), (|77|), ( f79| ) and (|90| ) we deduce 
Lemma 7.1 For u = PSi + P8 2 £ Z there holds 



(91) 



u(u x A g y +g x A u a ) = -8A d Id g(a) - 8A d R -ig(b) + O (e 2 | loge|) 



7.2 Mixed terms in P5\ and P5 2 

For u = PSi + PS 2 , we consider the first and the second integrals in (Q). We are interested in the terms 
involving both P8\ and PS 2 , namely 

f VPS 1 ■ VP8 2 + \ [ PSt ■ ({PS 2 ) X A (P5 2 ) y ) + ~ / P<5i • {{P5 2 ) x A {P8 X ) V ) 
Ja Ja J Ja 



+ - / PSr ■ {(PSi) x A (PS 2 ) y ) + - PS 2 - {{P5 2 ) x A (Ph) y ) 
•J Ja & Ja 

+ ~ l PS 2 - ((PSl) X A (P6 2 )y) + \ [ P6 2 - ((PSl) x A (PS^y) . 
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Integrating by parts it is easy to we see that the last expression becomes 

f VP6 1 • VP5 2 + 2 [ PS 2 - ((P5r) x A (PSi) y ) + 2 f PSf {{PS 2 ) X A (P6 2 ) y ) 
Jn Jn Jn 



in Jn 
Integrating by parts the first term we get 



/ VP5 1 ■ VP5 2 = - [ AP6 2 ■ P5i = -2 f P5 1 ■ ((S 2 ) x A ($,)„), 
Jn Jn Jn 

so we are left with 

2 / P8 2 {{P5 X ) X A {Ph) v ) + 2 I P6i ((Pda)x A (P5 2 ) y - (5 2 ) x A (5 2 ) y ) . 
J si Jn 

Lemma 7.2 For PS% + PS 2 e Z, there holds 

2 f PS 2 ((P8i) w A (PS^y) + 2 f PS! ((P5 2 ) x A (PS 2 ) y - (S 2 ) x A (S 2 ) y ) 



(92) 



2 / P* 2 ((*i)xA(*i) y ) + e(Ai,A a ). 
Jn 



Proof. Since the difference between the l.h.s. and the r.h.s. of 



2 / PSx ((PS 2 ) X A (PS 2 ) y ~ (S 2 ) x A (S 2 ) y ) + 2 PS 2 ((P8i) x A (PS t ) v - (Sx) x A (6^) , 
Jn Jn 

it is sufficient by symmetry to estimate one of the two terms in the last expression. We have 

Ph {{P5 2 ) x A (PS 2 ) y - (S 2 ) x A (S 2 ) y ) = J PSi ■ [(tp 2 ) x A {cp 2 ) y - {f 2 ) x A (5 2 ) y - (5 2 ) x A {tp 2 ) y } 
n Jn 

Using equations (|25|)-(|2^), choosing r < \C and setting Sl a ,b = fi \ (B T (a) U B T (b)), we find 



PSl ■ (((P2)x A (lfi2)y) 



PSi ■ {(cp 2 ) x A (6 2 ) v ) 



< 



< 



< 



C 
AiAo ' 



C_ 

C_ 

A 2 



\Ph\\{5 2 )y\+ \Ph\\{h) y \+ l^llK^I 

B T (a) JB-r(b) Jsi a , b 

C , c _logA 2 , C 



A1A2 A1A2 A1A2 



and an analogous estimate for the term P5\ ■ {(S 2 ) x A ((p 2 ) y ). This concludes the proof. 
Lemma 7.3 For P5\ + PS 2 S Z, there holds 



2 PS 2 - ({5i) x A 
Jn 



16 A fa'i-a'i dhi 
-r\i 



(93) 



A1A2 
16^o 



16An 



' - 22 



""1 _L dfl2 I h\\ 



' a x a 2 dhi \ \QA n ( ' a x a 2 dh 2 



+ e(Ai,A 2 ). 

PROOF. The left-hand side of ( |93| ) is given explicitly by 
(94) 

n 3 n 3 p 

2 PS 2 - ((6i) x A (6i) y ) = 2 ]T ry / (<5 2 ), • ((Si) x A (ft),,), - 2 £ r« / ■ ((Si) x A (fc),,)*, 

./fi j _-_ n -/n „• -/fi 



2(» 



where {rjA are the entries of the matrix R. We are now going to estimate these integrals. We recall 
that, by @ 

^(0 = (J^Mt, b ) + 0(A 2 - 2 ), 6 ) + 0(A 2 - 2 ), 1 - A/^, b) + 0(A 2 - 3 ) 

Taking this into account, we find that the terms in ( p4| ) involving the coefficients rn, ri 2 , ?"i3,r3i and 
are given respectively by 



(95) - 32A?A 2 

(96) - ;;2\yA 2 

(97) 16A 2 A 2 



(a; - b 1 )(x - ai) 



„ (i + A?ie - «i 2 ) 3 (i + Aiie - &| 2 ) 

:; N / (x-b 1 )(y-a 2 ) 

(l + A 2 |£-a| 2 )3(l + A 2 |£-6| 2 ) 



, ^,3 f (A 2 - 1 / tl (e,6) + Q(A 2 - 2 ))(x 1 -a 1 ) 
+ 1 L (l + A?|C-a| 2 ) 3 ' 



32Af 



3 /" (A 2 - 1 /i 1 (C,b) + Q(A 2 - 2 ))(y- a2 ) - 



(1 + AfK-al 2 ) 



2\3 



(^r- 6x) (1- A 2 |e-a| 2 ) 



(l + A 2 |e-aP)3(l + Ai|e-6| 2 ) 



16A 2 



(A 2 - 1 fa 1 (e,b)+Q(A 2 - 2 ))(l-A 2 |e-a| 2 ) 
( 1 + A 2| e _ a | 2)3 



(98) 
(99) 



16A? 
-8A 2 



n v l + A 2 |£-6| 2 A 2 



-^h 3 (^b) + 0(X^) 



(x - ai) 



(l + A 2 |e-a| 2 ) 3 ' 



1 + A 2 |e-6| 2 A 2 



-4^6) + 0(A 2 ^ a-^-P) 



(1 + A 2 |e~a| 2 ) 3 ' 



The terms involving the other coefficients of the matrix R can be estimated using the above ones, and 
will be taken into account later. 



Estimate of (|95|). Using the change of variables Ai(£ — a) i— > £, equation ( |95| ) becomes 

.A 2 /■ s(x-Ai<7i) , on /" (A2 1 /ii(Af 1 C + a,6) + 0(A" 2 ))x 



(100) 



-32- 



32 



A i K, (i + iei 2 ) 3 (i + fie - ah 2 ) ' ""V.. (i + l£l 2 ) 3 

We estimate the first term in ( |l00[ ). Consider the following subsets of the domain of integration 



We can write 



(101) 



i+lie-A^ 2 ) = (1 + AIM 2 )" 1 | J 



(1 + A||a| 2 ) 



On the set i3i we have the following inequality 



Ai 
V 



< 



A 2 |a| 2 , A 2 ^' 2 



16 



< -MM' 



hence, from a Taylor expansion, we obtain the following uniform estimate 



(102) 



f|£| 2 -2^-<7 



(1 + Aikl 2 ) 



- 1 



(1 + A 2 H 2 ) 



< C 



( 



ICI 5 



VAfkl 4 A 2 |^ 
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Using equation (1102) and some elementary computations we find 



x(x — AlCTl) 



(i + iei 2 ) 3 (i + f ie-AH- 



x(x — \\0\) 



(103) 

On the set £? 2 we have 

and hence we deduce easily 
(104) 



+ e(Ai,A 2 ) = 



(1 + 1^)3(1 + ^^2) 

A 1 



1 + AlkP 



A 



2al 



2 (l + Al^l 2 ) 2 (1 + A!M 2 )K 



+ e(Ai,A 2 ) 



T AfFF I 1 " 2 M 5 J + e(Al ' A2) - T^fH^ 



|?|>A 1 H-|e-A 1( x|> i A 1 | ( T| 



e(Ai, A 2 



\x - AiCTil 



< 



C 



a 2 (i + iei 2 ) 3 (i + f ie - AKrp) A >i 4 ^x^i i + f ki 



In B 3 we have |£| > and |£ — Ai<r| > ^+j^, and hence 



< C 



logAi 
Afkl 4 ' 



(105) 



[e||as — Aio-i| 



« 3 (i+iei 2 ) 3 



- C|fiil,8l WAS " C AfR 



Let us now treat the second term in (10C). Reasoning as above we find 



(106) 
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(A 2 - 1 /i 1 (A^ 1 e + a,&) + 0(A 2 - 2 ))a; 1 A 9ft 



+ 



2\3 



AiA 2 2 <9x 



-(a, 6) + e(Ai, A 2 ). 



Hence, using formulas (103)- (106), we are able to estimate (|95|), and we find 

16A / cr\ — erf dhi 



(107) 



2 / (P5 2 ) x A 



Ai A 2 



^(o,6)) +e(Ai,A 2 ). 



Estimate of (|96|). The proofs of the estimates of this and the remaining terms will only be sketched, 
since they are similar to that of (p5[). Using the usual change of variables, equation (p6|) becomes 



(108) 



-32^ 
Ai 



x(y - Aicr 2 ) 



i.. (i + iei 2 ) 3 (i + f ic-Aiap 
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(Aj 1 /n(Ar 1 e + Q,fe) + Q(A 2 - 2 ))y 
(1 + ICI 2 ) 3 



To treat the first integral in ( 108 ) we begin by dividing again fl^.a into the above sets B\, 2? 2 , B3. 
Reasoning as before and neglecting the higher-order terms we find 



A 2 



x(y - Ai<t 2 ) 



Ai J Bl (l + |eP)3 (i + ^- Xl a\ 

(A 2 - 1 /i 1 (Ar 1 e + a,&) + 0(A 2 2 ))y 
(1 + iei 2 ) 3 

Hence, using the last two equations we deduce 



<7lC 2 

~ A ° \ \ 1 14 +e(Ai,A 2 ; 
AiA 2 |cr| 4 

1 A dhi 
A^Tlfy 



(a, 6) + e(Ai,A 2 ). 



(109) 



2 / (P5 2 ) x ((Ji) x A (5i) tf ) 2 



AiA 2 



32^ + 16^(a,6)) +riA l .A,). 



2cS 



Estimate of (p7|). Wc turn now to the term involving the coefficient ri3, (|9^), which can be written as 
(1 - \tf)( Xl - Ann) A (^(A^ + a^J + OCAa^a-lSl 3 ) 



(110) 16^ 



A i -/n Al(0 (i + |e|2)3 ( 1 + || C _ Al(T p 



- 16 



Using (|27|), (102) and reasoning as in (103) one finds 



i + l?l 2 ) 3 U + ?ie-AH 



e(Ai, A 2 



Similar estimates hold if one integrates on the sets 2 and B3. Moreover, using (|27|) and elementary 
computations one finds 



(A^MA^ + a,b)+ Q(A 2 - 2 ))(1 ~ |£| 2 ) _ 
(l + lCl 2 ) 3 (l + f le-Axal 2 ) 

From the last two equations we deduce 

(HI) 2 / (PJ a ) x ((50* A (6x) y )a = e(Ai, A 2 ). 



e(Ai, A 2 ). 



Estimate of (|98|). The expression in (98) becomes 
32 



/ I i^(Ar^ + a, b) + 0(A- 3 ) * 



Reasoning as above, we obtain 

2 



^ 3 (Ar X £ + a, b) + 0(A 2 - 3 ) ] (1+ ^ |2)3 = e ( A i> A 2 



and that the integrals on the sets £? 2 and B3 are also of order e(Ai, A 2 ). Hence we find 
(112) 2 / (PS 2 ) 3 A = e(A 1; A 2 ). 



Estimate of (J99|) - We turn now to the term involving r 33 . Using the above change of variables, (99) 
becomes 

-4^ 3 (Ar 1 e + a,&) + o(A 2 - 3 ) ' 1 ^ 



■16 



1 



2 A9 



Using equation (|2^) and reasoning as above one finds 

(113) 2 / (P<5 2 ) 3 ((Si) x A (50^3 = e(Ai, A 2 ). 



(1 + ICI 2 ) 3 ' 



Other estimates. From the estimates of the terms (|9Jj)-(|39|) one can deduce also those involving the 
coefficients r 22 ,ri 2 ,r 2 3 and r3 2 . In fact, it is sufficient to permute the coordinates x and y in a suitable 
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way. Thus one finds 

(H4) 2 J^(PS 2h ((*). A (S 1 ) v ) 2 = + ^(a, 6)) + e(A 1; A 2 ). 

(H5) 2 JjPS 2 ) 2 A (50,)! = ^ + ^(o, 6)) + e(A 1 , A 2 ). 



(116) 



(P5 2 ) 3 ((5i) x A(8 1 ) y ) 2 = e(\ 1 ,\ 2 ). 



2 / (P*a)2((*i)«A(*i)»)s = e(Ai,A 2 ). 



Hence the conclusion follows from pfl), p§, (111), (112), p3), and pj-pfl) 



7.3 Expansion for k bubbles 



In this subsection we consider the case of k masses. When k = 2, from Proposition 5.1, and Lemmas 
(O) O, W3 we find 



16 



(117) 



+ 



9 

16A 
A1A2 



1 



1 r ~ r 



I e (u) = —A a + 8A ( ^Hia) + -^H(b) - —d Id g{a) - —d R -ig(a) 



Ai 



of - of dhi . ,A / of _ °1 



A 2 
2 



01(72 . \ / 01 0-2 . 9^2, ,v 



+ e(e, Ai) + e(e, A 2 ) + e(Ai, A 2 ). 



where u = P<5i + PS 2 and are the entries of the matrix R. 



We consider now the more general case of k masses. Given two bubbles 5, = Riir and 5k = Rk^ (we 
recall the definition of the stereographic projection 7r in Section ||), where Ri, Rk £ 50(3), we denote by 
Rik the matrix RT o R k . By invariance under rotation, it is clear that the interaction between Si and 8k 
is the same as the interaction between ir and R~ 1 8k = ilfcTr. 

In the expansion of the Euler functional for k masses, since I e is cubic in u, we are going to find mixed 
terms of the form J n PSi ■ (P5j A PSk), where i,j and k are all different. Since we are assuming that the 
distance of the points Pi,Pj and pk is uniformly bounded from below, there holds 



(lit 



[ PSi ■ {PS ) A PS k ) < C-^—, i^j^k,i^k. 



It follows that the interaction among three distinct bubbles in Z is negligible with respect to the inter- 
actions with g, SI and the interaction between two bubbles. 



We recall the definition of the quantity e(e, Ai, . . . , At) in Section g. Using equation (118), and 
omitting some straightforward but tedious computations we obtain the following Proposition. 



Proposition 7.4 Let C > 0, let k G N and let Z be defined by (plj). For i 7^ j let us set 



(119)F n (Pi,p j ,R i ,R j ) = mA 



(Rij)u 



(Pj -Pi)l - (Pj ~Pi)i . dhi 



\pj - Pi 



{Pj -Pi)l - (Pj -Pi)l . dh 



\P]-Pl\ l 



r r> \ 1 q (Pj-P 1 )i(Pj -Pih dh x \ ( Jpj ~p t )l(Pj -Pt) 2 dh 2 



and 



Sn, ff (e,pi, . . . ,Pk, Ai, . . . , A fc ,i?i, . . .,R k ) = ^Fn, g (pi, Xi,Ri) + V" Fn ^ £ ' Pl,Pj ' Rl ' R 

i < — ' AiAj 



i=l 



l<j 
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Then there holds 



I e (u) = -g-^o + ^n, g + e(e, Ai, . . . , A fe ); 



and 



whe 



dl s (u) = dZ n , g 
dpi dpi 



+ e(e, Ai,...,A fc ); 
dI E (u) 9£ 



dl £ (u) 9E n , s , 1 , . 



9A 4 



9A 4 Aj 



9ft 



^ + e(e, Ai,...,A fc ), 



9il 



i Afe); 



7.4 Some remarks 

In this subsection we consider the expansion for 2 masses with zero boundary data. Our goal is to 
extremize the functional in (117) with respect to a, b, ^ and R. Letting G denote the Green's function 
of f2 and setting 



dG 

Gi(o,0 = ^-(a,0; 



G 2 (a,0 = — (a,0 
da 2 



there holds 



<rf - of dhi , , , 



(7 



(71(72 . 9/l 2 , 

2 H^ + 9a>< 5) 



dx 
dG 2 



{a, by, 



dx 



(a, by, 



<Jio- 2 dhi 

2 TIT + ^r~( a ' 6 ) 

|cr| 4 dy 
of - a\ dh 2 



dy 



dG x 
dy 
dG 2 



(a, b); 



+ —±(a,b)\ = —±(a,b). 



dy 



Using these expressions, the expansion of I s (u), with u = PSi + RP8 2 becomes 



16 



9 

16A 
A1A2 



1 / dh\ dh 2 



I e (u) = -A + 8A _ _i + _± Ua,a) + I -± + -± (b,b) 



X\ \ dx dy 



1 / dh\ dh 2 



\\ \ dx dy 



1 i\ . dGi dG 2 dG 2 

rn-^(a,b) + r 12 -^-(a,b)+r 21 — (a,b)+r 22 —(a,b) 



+ e(Ai,A 2 ). 



The entries 7^3 and of the matrix R appear as lower order in the above formula, see Remark |5.3| (b) . 
We can write 

dG dG dG dG 

m— l(a, b) + ri 2 —±(a, b) + r 2 i—^(a, b) + r 22 ^(a, b) = ei ■ flVGi(o, b) + e 2 • i?VG 2 (a, b). 
dx dy dx dy 

As in [pl[ , Lemma 5.4, the extremization with respect to R gives 

ei • RWGi(a, b) + e 2 • i?VG 2 (a, 6) = ± (|VGi(o, 6)| 2 + |VGi(a, 6)| 2 ± 2|VGi(a, 6) A VG 2 (a, 6)|)* . 
Hence, setting iT = (hx) x + (h 2 ) y , we are left with 
16 



h{u) 



± 



9 

16^o 
A1A2 



1 ~ 



A + 8A [ -^H(a) + ^H(b) 



1" 



ggi 9G 2 
9a; 9y 



(a, b) + 



( dG 2 dd 
\ 9a; 9y 



(a, 6) 



e(Ai, A 2 ), 



where the + and — signs inside the square brackets are opposite (hence there are four different pos- 
sibilities). I e (u) has the form c + an£ 2 + a 22 ^ i 2012^1^2, with > 0. Thus if we consider the 
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case c + an£i + 022^2 + 2ai2£if;2, we notice that minimizing X) a ijCiCj'/ICI we necessarily need to se- 
lect £ = (£1,62) with £1^2 < and so this case does not arise. Thus the only case that remains after 
cxtrcmizing is 



16 



(120) 



16^o 
A1A2 



dG 1 dG 2 
dx dy 



(a, &) 



( dG 2 dGi 
\ dx ~^ dy 



(o,6) 



+ e(Ai,A 2 ). 



8 Proof of Theorem \T72\ 



In this Section we prove Theorem [l.2[ We begin with the following Lemma, proved in |26 and which 
follows from straightforward computations. 

Lemma 8.1 Let ui 6 (0, 1), and let a w = (ui, 0) € D. Define also : <9£> — » K 3 as 

a; — w y 



r,0 



> (a; — w) 2 + y 2 ' (x — to) 2 + y 2 
Then, letting g^ be the harmonic extension on D ofg^, there holds 



(121) 

and 
(122) 



9w{x,y) 



x — lu(x 2 + y 2 ) 
(1 - ujx) 2 + (ujy) 2 ' (1 - ujx) 2 + (coy) 2 



,0 



|V^| 2 + 2|(^) ;c A(^) J/ | 



((1 -w.t) 2 + uj 2 y 2 ) 2 ' 



(x,y)€dD. 



[x,y) e D, 



(x,y) e D. 



By equation (66), the concentration points of blowing-up solutions as e — > are critical points of the 



function Vg l ±2 lg* A g« _ j n f-h e nex ^ Lemma we describe the critical points of this function in the case of 



Lemma 8.2 Let ui € (0, 1), and let g u be as in Lemma 8.1. Then one has 



:= 



(\Vg\ 2 ±2\g x Ng v \ 



= V2-, 



(1 - x 2 - y 2 ) 



V H J ' "(1 -lux) 2 + {ujy) 2 ' 

The point (uj, 0) is a non-degenerate global maximum for W u and 

(123) W u (u,0) x2 
The Hessian of W u at (lo,0) is given by 

(124) D 2 W u (u),0) = -2V2 ( TT 



1 -uj 2 



(1-u 



Remark 8.3 From equation (124), the fact thatV g(uj,0) 7^ 0, and from Theorem D in \2Ql it follows that 
problem (^[) admits a solution concentrating at (lu,0) as e — > 0. The image of these solutions converges 
to a sphere of radius 1 ce ntered at ( 0,0,-1), since <p a ,x — * (0,0,-1) as e — > 0, see (f22|). 

Note that, from ( 123 ) and ( 124 ), attains a sharp maximum with highly non- degenerate hessian 
when uj is close to 1. We will use this fact to glue k single bubbles showing that, for a suitable boundary 
datum, the interaction of this datum with the bubbles is stronger than the interaction among different 
bubbles. In the next Lemma we give quantitative estimates of the gradient of ' Fo igu (see Proposition 5.1) 
in a suitable neighborhood of one of its critical points. 
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It is classical to represent a rotation Rq £ 50(3) using the Euler angles in the following way 
Ro = 



cos ip cos 4> — cos 9 sin sin ip cos ip sin + cos 9 cos sin ip sin sin I 
- sin ^ cos <p — cos 6* sin cos "0 — sin ip sin + cos 9 cos cos ■0 cos ip sin ( 
sin 9 sin — sin 9 cos cos 6 



where 9 G (0, 7r), 0, G (0, 27r). For us it is convenient to use coordinates different from the Euler angles, 
in order to have a smooth parametrization near the identity matrix. A rotation R will be parameterized 
as 



(125) R- 1 
namely as 



cos ip cos — cos 9 sin sin ip cos ip sin + cos cos sin i/j sin sin 9 

— sin 6> sin sin 6* cos (p — cos 6> 

— sin ij) cos — cos 9 sin cos ip — sin sin + cos cos cos cos sin ( 



R- 1 



'10 
-1 | R . 
,0 1 



i? is the identity matrix for 9 = ^, and cf> = ip = 0, and the angles 9,ip,(p are smooth coordinates near 
the identity. In fact there holds 



dR- 



09 



'0 
= I -1 
,0 1 



dR- 



dip 




dR- 



00 




when 9 = ^, = = 0. We will show that the identity matrix is critical with respect to the rotations 
for the quantity dj^-i g w (w, 0). There holds 

dn-iQoj = (cos ip cos 6 — cos 9 sin sin-0) -f (cos -0 sin + cos cos sin 0) ^"^ 2 

ax ox 

— (sin sin 0) — h (sm6> cosc^ 



dy 

From simple computations one finds 



Ox 

and hence 



(126) 



dy 



(1 — lux) 2 — Lu 2 y 2 



dy 
d{gu>)i 



(1 — ujx)usy 



((1 — lux) 2 + Lu 2 y 2 ) 2 ' ch/ 
(cos cos — cos sin sin -0 + sin 9 cos < 



cte; ((1 — lux) 2 + Lu 2 y 2 ) 2 ' 

(1 — lux) 2 — Lu 2 y 2 

2(cos ip sin + cos 9 cos sin ip + sin 6* sin < 



((1 - wi) 2 + u 2 y 2 ) 2 
. (1 — ux)ivy 



((1 - cjx) 2 + w 2 ?/ 2 ) 2 



In the next Lemma we study the critical points of Frj l9ul for £ ~ (w,0), i? ~ Id, A ~ 2e 1 and £ small. 
We use below the coordinates 9, ip, in (125) to parametrize the matrix R. 

Lemma 8.4 Let lu G (0, 1) and let g^ be as above. Then, for fixed e, the point x = lu , y = 0, X = -, 
9 = 77, tp = 0, = is critical for Fu, 9lo (e, £, A, 9, ip, 0). For \x > define the set 



\x-lu\ <^{1-LU 2 ),\ V \ < f l(l-LU 2 ), 



x- 2 - 










•-f 


< 


e 


e 







<fi,\ip\ </i,|0l <M 



Then for fi sufficiently small and lu sufficiently close to 1, there exists a universal constant Co independent 
of e, \x and lu such that 



(127) 



VF D ^(e, x )-x>C 



F 2 II 2 
1 £ f* 



on 970, and hence deg(V-FD )9u , 70, 0) = 1, 



(1 - cj J 

where \ denotes the set of variables x, y, X, 9, ip, 0, and the gradient is taken with respect to \- 
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Proof. We recall that the functional Fjj i9ui is defined by 



where 



(i-iei 2 ) 2 ' 



and where dR-ig w (£) is given by (126). In particular there holds 



2 7T 



e' 2 



VF D „ e,w,0,-,-,0,0 =0; HessFjj „ e,w,0, -,0,0 =2 



2 7T 



rA 



where 



/ 1 , 3uj 2 



A, 



" 71 — "TIT T 71 — 777\7 



V 



' 1 (1-w 2 ) 













1 euj 
'2 (1-w 2 ) 



±e 2 






2 '"'7 ~(l-w 2 ) 2 
\ 



0- 



* ? 

3 









5 / 



We point out that the matrix A u , £ is positive-definite and non-degenerate. Using simple but tedious 
computations, one finds 



(128) 



/j, 3 e 2 



\VF D , 9u , (e, X )-X- A,eX\ < C jfrjlf 



for x,y,X,6,ip,(f) € 



Then the conclusion follows from the fact that cj,0, |, f ,0,0 is a critical point of Fu, gio , from (125), and 
from the fact that e is positive definite. ■ 



Now we are in position to prove Theorem 1.2. For the main idea see Remark p.3| 



PROOF of Theorem |L2| Let A = {Si U • ■ ■ U S k } be as in Theorem [L2j and let {vi, . . . ,v fc } C K 3 
denote the centers of S%, . . . , Sk respectively. Note that, since all the spheres have radius 1 and since they 
all pass through the origin, one has |vj| = 1, for all j = 1, . . ., k. Let TZi, . . . ,lZk £ SO (3) satisfy 

(129) 7^(0,0,-1) = vj, for all j = l,...,k, 

see Remark 8.3. Let lu S (0, 1) and define Gk,u '■ dD — > R 3 by 



3=1 



(x,y)edD, 



where 



27Tj~ 

cos I x - 

k 



2-KJ 

sin — ] //. 



sm — 



It is clear that the harmonic extension Gfc jtlJ of Gfc jU to the interior of D is given by 



2ttA 
k ) 



Gk,u>(x,y) = ^Tlig jtU {x,y) 

3 = 1 

(130) = J2 n ^(( cos ^r- 



3 = 1 



2nj\ ( 2nj s 
sm— y,- sm— ,,- 



cos^)?/), {x,y)eD. 



where g^ is given by (121). Our goal is now to study the critical points of the functional YtD,G k „ defined 
in Proposition 7.4. 
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We introduce coordinates 9j, ipj and (pj parameterizing a rotation Rj (note that this is a generic 
rotation, which differs from the fixed rotation IZj) in the following way 



(131) 



cos ipj cos ipj — cos 9j sin (pj sin ipj cos ipj sin <pj + cos 9j cos 0j sin ipj sin -0-,- sin 9j 



cos f 



sin cos 4>j — cos #j sin (pj cos V'j — sin ipj sin + cos 9j cos </>j cos ipj cos i/^- sin 0j 



RfTlj. 



The choice of this parametrization will become clear below. Note that for Qj ~ ? and ipj,(pj close to 0, 
these angles are a smooth parametrization of SO (3) near 7?.j. Define also the set 



, 27TJ . 2%j 

[Xj ,yj)-u\ cos — , sin — 



< /Lt(l - LO 2 ), 



A, 



< 



/' 



< H, \ipj\ < |03 1 < M 



We are going to prove that for e sufficiently small and for ui sufficiently close to 1, the functional Ed 
has a critical point with Xj, yj, Xj, 9j, ipj, <pj G for all j = 1, . . . , k. 
By Proposition 7.4 and by the definition of FD,G k u we have 



(132) 



k ( ) 

3=1 i<j ' ^ 



*<3 



By invariance we can write 



A fe L fc 



H(£ k )-eXd R -iK k g u (t,) 



(133) 



= F D<gui (e,p k ,X k ,TZ k R k ) = F D ^{e,x k ,y k ,X k ,9 k ,ip k ,(p k ). 



We remark that the function in (133) is exactly the one studied in Lemma S.4. This justifies the choice 
of the coordinates 9j, ipj, <pj in (131). 
There holds 



d 



a — p D, 9a (e, x k , y k ,X k , 9 k ,ip k ,(p k ) - 8—A Q V" - — d R -ig^^[p k ) 
ox k \ k 'rr ox k k 



E 

l^k 



d F D (pi,p k ,Ri,R) 
dx k 



XiX 



Sincc all the A, 's are of order e 1 , and since the mutual distance between the points pj 's is bounded from 
below, it is easy to check that for p, sufficiently small 



(pi,Xi,Ri) G Ty i 



£ 



±-d R -,g^{p k ) 



d F D (pj,p k ,Rj,R k ) 



dx k 



XjX k 



<Ce 2 j^k, 



where C is a positive constant independent of e, lu and fi. Hence the last formula and (132) imply 



d 



■^D,G k 



d 



dx k ' dx k 
Using similar estimates we find 



FD,g w (e, Xk , yk , Afe , 9 k , ip k , <p k ) 



< Ce 2 if (p 3 , Xj , Rj) G 7i for all j = 1 , 



,k. 



(134) 



C^D,G k 



^(F D , gu ,(e,x k ,y k ,X k ,9 k ,ip k ,4> k )\ < Ce 2 ; 



afr s AG fc ,„ - -£-F D , gu ,(£,x k ,y k ,X k ,9 k ,ip k ,(f> k ) 



<Ce 3 , 
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provided (pj, Xj, Rj) 6 for all j = 1, . . . , k, Here £ denotes the set of variables Xk, TJk, Ok, ipk, <t>k> an d 
where C is a positive constant independent of e, u> and fi. 

Let us fix ^ and e sufficiently small such that (127) and (134) hold. Then we have 

V x £ AGf _ • X > Co' 1 (1 1^2)2 - Ce2 > * e if (Pi. A i> R j) e ^ fOT a11 3, 

where C and Co are independent of e, u> and /i. Now, choosing u> sufficiently close to 1, depending on C, 
Co and /x, and reasoning in the same way for the indexes different from k we obtain 



V X3 Xd,g>, u • Xj > -^- £ 2 , deg(V Xj E D>GfciB , 7-?, 0) = 1, if (pj,Xj,Rj) G 7^ for all j, 



2 

where denotes the set of variables Xj,yj, Xj, 9j, ipj, 4>j. For the above choices of \i and u>, let I Gfc ^oj.e 
denote the Euler functional I £ corresponding to the boundary datum Gk^(n)- By Proposition 4J3, for e 
sufficiently small we obtain 

C" 1 

V v ln.<;. .. • X 3 > ^-£ 2 , deg(V Xj / AG; _, e , T£, 0) = 1, if [pj^Rj) S T^f for all j. 



Then, by Proposition 4.3 and Lemma |8.5| below, letting e — > 0, we find a family of solutions u £jM of 
/ Gfc = satisfying, up to a subsequence 

u e ^(D) — > A M = {Si,^, . . . , S'fc.p} in the Hausdorff sense as e — > 0, 

where Si,^, . . . , Sfe 1(U are spheres of radius 1 passing through the origin and lying in a neighborhood of 
order /i of Si, . . . , Sk respectively. Now we can choose /z(e) — > sufficiently small as e — > 0, and find a 
corresponding w(e) — > 1 such that the solution ite !Al ( e ) of J Gj , ^ = obtained with the above method 

satisfies 

u e ^{D) — ► A = {Si, . . . , Sfc} in the Hausdorff sense as e — > 0. 
This concludes the proof of the Theorem. ■ 

Lemma 8.5 Let g : dfl — > R 3 6e a smooth function, let k € N, C > 0, and Ze£ Z 6e defined as in (|3 
-Lei u be a solution of (P s ) of the form 

k k 

u = PRiS Pu \ i + w; with PRi8 Pi ^\ i € Z, and ||u>||jji(f2) ^ as e — > 0. 

i=l i=l 

T/ien |MU°°(fi) -> as e — > 0. 

PROOF. In the following we simply write Si for Ri5 Pu \ i , and we let ^ be the function in ( [H| ) corre- 
sponding to 5j. The function w satisfies 



Aw = 2 -cpi) + w + eg) x A (E a -(<*i - Vj) + w + - Ei(*i)« A in fi > 

w = on <9f2. 



where g, as before, denotes the harmonic extension of g to fl. Expanding the wedge produce on the 
right-hand side we obtain (as before PSi = Si — (fi) 

Aw = 2 ^(P^)* A (P^)j, + [(P8%)* Affi 9 + Wl A (P^) y ] + 2e ^ [(PSi) x A g y + g x A {PSi) y } 

ijtj i i 

(135) - 2 [(dj) x A (<pi) y + (<pi) x A (Si) y ] + ^2((fi)x A (</?i) a + 2w x A w y + 2e(w x A g y + g x A w y ) 

i i 

+ £ 2 9xAg y . 



3G 



Figure 4: the points pj on the disk D and the bubble generated by g u j 

Using (p5|), ( ptf ) and some elementary computations, for any p > 1 the first term in the right hand side 
can be estimated in the following way 



IKiMiJxA^OvlU'Cii) <C(C,p) 



1 




From standard elliptic estimates it follows that 

\\(Ar 1 ((P6 l ) x A(P6 ] ) y )\\ L ^ m <C(C,p)ei, 
where (A) -1 denotes the Green's operator for A in SI with Dirichlet boundary conditions. Let us focus 



now on the second term in (135). Writing for brevity ip = PSi, one has 

(ip x A w y + w x A ipy) = [J(w 2 ,ip3) + J{ip2,w 3 )]i + [J^i,^) + J{i>i,wa)]j + {J{wi,^) + J(ipi,w 2 )]K 
where J(F, G) = F x G y — F y G x is the Jacobian function. By the result in |Q there holds 
(136) 11(A)" 1 (ii x A w y + w x A ip y )\\ Lao{sl) < C\\P5i\\ H i m \\w\\ H i {n) as e 0. 



The remaining terms in (135) can be estimated as in (|13G|) 



Remark 8.6 FFii/i an easy modification of the above arguments we can easily obtain the limit configura- 
tion {Si, . . . , Sk} with a boundary datum of the form eG, for some fixed function G on 3D independent 
of e. 



Remark 8.7 We remark that to obtain L°° estimates on the solutions of (\P e \) , we use in a crucial way 
that these solutions satisfy the H -surface equation with H = constant. Such estimates are not available 
for general Palais-Smale sequences, as exhibited in J7^/ . 

In Figure |i] we indicate the location of the points pj in D when w is close to 1, see the definition of 
7^. We al so p lot the boundary datum egw,j, which lies in a plane, and the corresponding bubble (as 
in Remark |8.3| ) whose center Vj which is, roughly, perpendicular to the plane of g u ,j- We note that the 
image of g^ is a great circle (the Kelvin inversion of dD w.r.t. the point (u>, 0)). In Figure |E| we plot the 
configuration of bubbles generated by the function eGk,u- Each bubble is nearly perpendicular to some 
gj,uj (whose sum is Gt-, w ). 
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Figure 5: the boundary datum eG k ,u and the corresponding configuration of spheres 

9 Appendix 

This Appendix is devoted to the characterization of the solutions of the equation / (S a ,x)[w] — 0. We 
can suppose by invariance that a = and A = 1, and we set <5 a ,A = 7T, see Section [|. If a function w 6 D 
satisfies l"(S)[w] = 0, then it solves the linearization of (|l5|), namely 

(137) Aw = 2 (w x AS y +8 x A w y ) , in M 2 , w S X>. 

After inverse stereographic projection, equation (137) can be equivalently viewed on S* 2 as follows 



(138) A go w = 2(sinipy 1 {we A 5^ + 5 e Aw v ), in S 2 , w £ H 1 (S 2 ;M 3 ). 

where (9, ip), < 9 < 2ir, < ip < ir, are spherical coordinates on S 2 and A go is the Laplacian with 
respect to standard metric on S 2 . 



To analyze (138) we shall use some properties of spherical harmonics that we now recall. Let P n {x), 



x e (—1, 1), denote the n-th Legendre function. We define the associated Legendre function Pn{x) by 

(139) pk {x) = { _ 1) k {1 _ x 2 ) %jjL_p n{x)) fc > _ 
The spherical harmonics are defined by 

(140) Y n , k (e,(p)=c nt \ k \Pi k \(coB<p)e M> , - n <k<n, 
where the normalization constant c nj |fc|, see equation (21), p. 171, is given by 



(141) Cn t \ k \ 



2n + lV /(n-|fc|)! 
Air J V (n + |fc|)!' 
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From ( |139| ) we easily have, by differentiation 

(142) P k+1 (cosip) = - sin tp(P£)' '(cos ip) - m cot tpP% (cos ip), 
and from equation (41), p. 107 in 

(143) P k+2 {co&Lp) + 2(m + 1) cot (cos ^ + (n - fc)(n + fc + l)P^(cosy>) = 0. 



From (141) we also have, for 0<k<n,n>2 



(144) 



C n ,k+l 



< v3/2n, e n>k ■= (n — k)(n + k + 1) 



Cn,fc 
Cn ,fc_l 



In the sequel we will often write P„ (<^>) for P* (cos <p) . 

We define the finite-dimensional subspace Ti n to be the linear span of Y n<k , —n < k < n, and Ti n the 
subspace of L 2 (S 2 ; K 3 ) consisting of vectors it; = {fx,f2,h) with e P„, i = 1,2,3. Recall that any 
function w £ L 2 (S 2 ; K 3 ) can be decomposed orthonormally asw = S^Lo w ™' w ^ n w « € Ti n - We have 



Lemma 9.1 Let w : S 12 

io„ = for n > 4. 



&e a solution of (138), and let w = X)^Lo u '™' with w n £ 7i„. TTiera 



Proof. We first claim that for any F £ Ti n there holds 



(145) 



r(F) :=A 90 P- 



sm ip 



(F e A S v + 5g A FJ e 7U- 



If our claim is verified, to prove the Lemma it will be enough to pick a solution w to (138) in Ti n and to 
show that w n = for n > 4. 

We now prove our claim. W.l.o.g. pick F £ Tin of the form 



F = (a fcl c„, fcl P^(^) e ^fcl ^/? fc2 c Il , fc2 P^(^)e ^fc2 ^7fc3Cr ! , fc 3^ 3 (^)e ^fc3e ) = (a kl Y nM , f3 k2 Y nM ,~f k3 Y nM ) 

Next we have 6 V — (cos ip cos 8, cosy? sin 0, — simp) and Sg = (— sin ipsm9 1 sunp cos 9, 0). We will show 
that T(P) = -n(n + l)F - 2v, where 



— ik 2 Pk 2 YnM + l/2dn,fc 3 7fc3^ri,fc 3 + l — l/2e ni fc 3 7fc 3 5 / n,fc 3 -i 

ifeiafci5^»,fci - «/2(i„,fe 3 7fc 3 l / n^ 3+ i - i/2e n ^ 3 jk 3 Y n: k 3 -i 
i/2d n ^k 2 (3k 2 Yn,k 2 +i — i/2 e n,fcjAj^n,fc 3 -i — l/2d n! fe 1 afc i y„ ! fe 1+ i + l/2e„ i fe 1 afc 1 Y n ^ 1 _i 



(146) v 

Since v £ Ti n , our claim follows. It is evident that A„ P 



-n(n + 1)P, thus it is enough to show that 



the second expression in (145) is 2v. This follows by noting that 



-ik 2 /3k 2 Cn,k 2 sin tpP k2 (tp)e lk20 - ik 3 g k3 c nM cos ip sin 6 P* 3 (ip)e ik3 ^ 
ik 3 gk 3 Cn,k 3 cosipsm6P% 3 (ip)e lk3e + ikia kl c nM sinipP kl (ip)e ikl6 
ikia kl c nikl cos cp sin 9 P^(tp)e ikl0 - ik 2 /3 k2 c nM sin pP k2 (ip)e lk2f> 

-gk 3 c n ,k 3 sm 2 ipcos9(P k3 Y(p)e^ e 
~lk 3 c nM sin 2 ipsin9(P k *y(<p)e ik ° e 
K (3 k2 c n .k 2 sin 2 psmO(P k2 )'{ V y k29 + a kl c n , kl sin 2 ^cos9{Pt)'(cp)e^ 



FgA5 v 



8g A Ft* = 



Using ( |142j ) and ( pL43|) it is easily verified that (sin^) -1 [Fg A 6 V + Sg A F v ] = v. Then from ( |146[ ) and 
integrating the equation T(w) ■ w = on S 2 we find 



-ikf3 k a k + l/2d„ ifc 7fcO! fc+ i - l/2e„ ! fc7fe_ Q; fc-i 
= -2^( ika k (3 k - i/2d n , k j k f3 k+1 - i/ 2e n,*7*Pfc-i 

v i/2d n ,fc/3fe7fc+i - i/ 2e n,kPk7k-i - l/2dn,fcafc7fc+i + V 2e n,fcQ!fe7fc-i 
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where w = J2k( a k Y n,k, PkY n ,k-,lkYn,k) and where A 2 = J2\ a k\ 2 , B 2 = J2\Pk\ 2 , C 2 = J2 \lk\ 2 - Using 



(144), (147) and the Cauchy-Schwartz inequality we find 

i(n+ l)\(A 2 ,B 2 ,C 2 )\ < 2n (AB + y/3/2AC, AB 



< 2n 



V6 



r Zl2BC^[Zl 7 i{BC + AC) 



(A 4 + B 4 + C 



4x1/2 



Thus 



1 < (18 + V24) 2 , which implies n < 3. 



Lemma 9.2 The solutions of equation (138) are of the form 



ax 2 + (3x3 



'X\ 



w 



c + —axi + 7x3 + (a'xi + f3'x 2 + j'%3 ) x 2 



-f3xi - 7x2, 



K x 3, 



where c € R 3 and a, j3, 7, a! , /?', 7' € K are arbitrary. 



Proof. We denote by J(w) the r.h.s. of (138). From the proof of Lemma 9.1 it follows that J preserves 
the degree of spherical harmonic functions. Equivalcntly, J preserves the degree of polynomial functions 
in K 3 restricted to S 2 . By this reason and by Lemma 9.1, we can confine ourselves to study J just on 
polynomials of order 1,2 and 3. Since the computations involved in the proof are straightforward, we 
just give a simple sketch below, omitting some details. 
Using simple computations, we obtain 



(148) J(x u 0,0) = (0,2x2,2x3); J(x 2 , 0, 0) = (0, -2xi, 0); 
With a permutation of coordinates one also finds 

(149) J(0,x 2 ,0) = (2x l7 0,2x 3 ); J(0,x 3 ,0) = (0,0,-2x 2 ); 



■7(13,0,0) = (0,0,-2^). 
F{0, Xl ,0) = (-2x2,0,0). 



(150) F(0,0,x 3 ) = (2xx,2x 2 ,0); F(0,0, Xl ) = (-2x 3 ,0,0); F(0,0,x 2 ) = (0,-2x 3 ,0). 
Hence, letting w = (aiXi + a 2 x 2 + a 3 x 3 , 61X1 + 62X2 + 63X3, c\X\ + c 2 x 2 + c 3 x 3 ) we find 

(-61X2 + &2X1 - C1X3 + c 3 xA / aiXi + a 2 x 2 + a 3 x 3 \ 

aix 2 - a 2 xi - c 2 x 3 + C3X2 ; Aw = -2 61X1 + 6 2 x 2 + 6 3 x 3 , 
aix 3 - a 3 xi + 6 2 x 3 - 63X2 / \C1X1 + c 2 x 2 + c 3 x 3 / 

The system of equations Aw = J(w) admits the following solutions 
(151) 

with a, (3, 7 arbitrary real numbers. 

Let us now consider the homogeneous secon d or d er po lynomials. We have Axf = 2(1 — 3xf) and 
A(xiXj) — —GxiXj. Using the Leibnitz rule and (14S)-(150), we can compute J(w) when w has the form 



/ dl 


a 2 


03 \ 




b 


62 


63 






c 2 


c 3 ) 






w 



1 ' a\x\ + a 2 x\x 2 + a 3 x 2 + 04X1X3 + a 5 x 2 x 3 + a 6 x§\ 
b\x\ + 62X1X2 + 63X3 + 64X1X3 + 65X2X3 + 6 6 x| 

C\x\ + C2X1X2 + C3X2 + C4X1X3 + C5X2X3 + C 6 X§ J 



From the relation Aw = J(w), and using elementary computations we obtain 

w = (ax\ + [3x 2 + 7x3) + (S, T], a), 
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where a, (3, 7, S, rj, a are arbitrary real numbers. 

Let us now turn to the third order polynomials. We have 



Ax i — 6xi(l — 2x1); ^(xfxj) — 2xj(l — 6x 2 ); A(xiXjXk) — — !2xiXjXk- 



Again, the values of F on the third order polynomials can be computed with the Leibnitz rule and 
(p48|)-(^50|). Letting 



' a\x\ + a 2 x\x 2 



w = b\x{ + b 2 x{x 2 



c\x\ 



C 2 X\X 2 



a 3 x\x 3 - 
I- b 3 xjx 3 
h c 3 x\x 3 



a±x\x 2 x 3 

- b±x\x 2 x 3 

- C±X\X 2 X 3 



+ a$x\x\ + a§x\x\ + ajx^ 
+ 65.T1.T2 + b§x\x\ + b-jx 2 - 

+ C5X1X2 + CqXiX 3 + C-;x\ - 



- a§x\ + agx 2 x 3 + 010^2^3 ' 
b s xl + bgx%x 3 + bi x 2 xl 
c 8 xl + c§x\x 3 + c w x 2 x 2 3 



Hence, equating the coefficients in the above two expressions we find the following system and considering 
the equation Aw — J(w), we find a system decoupled in four parts. The first part consists of seven 
equations involving the seven terms a 2 , 0,7, aio, 61,65,66 and C4. Using simple computations one finds 
a 2 = 07 = aio = a, 61 = 65 = 6g = —a, C4 = 0, for some a € M. 

The second part consists of seven equations involving the seven terms a 3 ,ag,ag, 64 and c\,c^,cq. 
Using simple computations one finds a 3 = ag = ag = f3, ci = C5 = cq — —j3, 64 = 0, for some (3 £ R. 

The third part consists of seven equations involving the seven terms 04, 63, 6s, 69, and c 2 ,cj, cio- Using 
simple computations one finds 63 = 6s = 69 = 7, c 2 = c-j = Ciq = —7, 04 = 0, for some 7 g K. 

The fourth part consists in nine equations involving the terms a\, as, a§, b 2 ,b^, 610 and c 3 , eg, eg. Using 
simple computations one finds ai = 05 = ae = b 2 = 67 = 610 = c 3 — eg = eg = 0. 

The solution obtained in this way represent just the linear functions in (151), taking into account of 



the identity x\ 



= 1 on S . This concludes the proof. 



Proof of Proposition 4.1. Coming back to the space V, and using some elementary computation, 



the proof of the last statement follows immediately from Lemma 9.2. The first inequality is immediate 
to check. The second inequality follows from the proof of Lemma 9.2 when v S (B n >4Hr 



When v has 



some non-zero components in ©noHn, then it is sufficient to use straightforward computations, since 
we have to deal with finite combinations of spherical harmonics. Alternatively note that, since S is a 
mountain-pass critical point of /, the linearized operator possesses only one negative eigenvalue (with 
corresponding eigenvector <5), hence if v _L 6 and v _L Kerl (6), v must be a combination of positive 
eigenvectors. ■ 
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